The Astrophysical Journal 

Preprint typeset using KTgX style emulateapj v. 1 0/09/06 



THE M-a AND M-L RELATIONS IN GALACTIC BULGES, AND DETERMINATIONS OF THEIR INTRINSIC SCATTER 

Kayhan Gultekin', Douglas O. Richstone', Karl Gebhardt^, Tod R. Lauer^, Scott Tremaine"*, M. C. Aller^, Rale 
Bender^, Alan Dressler^, S. M. Faber^, Alexei V. Filippenko^, Richard Green'", Luis C. Ho^, John Kormendy^, John 

Magorrian", Jason Pinkney'^, and Christos Siopis'^ 

Received 2008 December 27; accepted 2009 March 26 

ABSTRACT 

We derive improved versions of the relations between supermassive black hole mass (Mbh) and host-galaxy 
bulge velocity dispersion (o) and luminosity (L) (the M-a and M-L relations), based on 49 Mbh measurements 
and 19 upper limits. Particular attention is paid to recovery of the intrinsic scatter (£o) in both relations. We 
find log(MBH/M,. ) = a+ plog(o/200 kms"') with (a, p, eo) = (8.12±0.08, 4.24±0.41, 0.44±0.06) for 
all galaxies and (a, p, Eo) = (8.23 ±0.08, 3.96 ±0.42, 0.31 ±0.06) for ellipticals. The results for ellipticals are 
consistent with previous studies, but the intrinsic scatter recovered for spirals is significantly larger. The scatter 
inferred reinforces the need for its consideration when calculating local black hole mass function based on the 
M-a relation, and further implies that there may be substantial selection bias in studies of the evolution of the 
M-a relation. We estimate the M-L relationship as log(MBH/M0) = a+ piog(Lv/10"Lj,v') of (a, p, Eo) = 
(8.95±0.11, 1.11 ±0.18, 0.38±0.09); using only early-type galaxies. These results appear to be insensitive 
to a wide range of assumptions about the measurement errors and the distribution of intrinsic scatter. We show 
that culling the sample according to the resolution of the black hole's sphere of influence biases the relations to 
larger mean masses, larger slopes, and incorrect intrinsic residuals. 

Subject headings: black hole physics — galaxies: general — galaxies: nuclei — galaxies: statistics — stellar 
dynamics 



1. INTRODUCTION TO BLACK HOLE MASS RELATIONS 

Studies of elliptical galaxies and spiral bulges ("hot" galax- 
ies) have led to the discovery that most such galaxies contain 
massive dark objects at their centers, presumably black holes 
(BHs; Kormendy & Richstone 1995; Richstone et al. 1998; 
Kormendy & Gebhardt 2001). Moreover there is a remarkably 
tight correlation between the BH mass and the slit-averaged 
velocity dispersion of the hot component of the galaxy (Fer- 
rarese & Merritt 2000; Gebhardt et al. 2000a). This M-a re- 
lation suggests a strong link between BH formation, galaxy 
formation, and active galactic nuclei (AGNs). 

The slope of the M-a relation has been estimated sev- 
eral times in the last ten years: originally 3.75 ±0.3 (Geb- 
hai-dt et al. 2000b) and 4.8 ±0.5 (Ferrarese & MeiTitt 2000), 
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then 4.72 ± 0.36 (Memtt & FeiTai-ese 2001b) and 4.02 ± 0.32 
(Tremaine et al. 2002), and more recently 4.86 ± 0.43 (Fer- 
rarese & Ford 2005) and 3.68 ±0.42 (for barless galaxies; 
Graham 2008). The reasons for the differences among the dif- 
ferent measures are discussed by Tremaine et al. (2002) and 
Ferrarese & Ford (2005). 

There is also a relation between the BH mass and the bulge 
or spheroid luminosity of the galaxy (e.g.. Dressier 1989; Ko- 
rmendy 1993; Magorrian et al. 1998, and the ratio of BH mass 
to bulge mass was found to be 2.2^q g x 10^^ by Kormendy & 
Richstone 1995), but the scatter in the relation is larger than 
in M-a. No other single parameter or combination of pa- 
rameters of the host galaxy has been found to predict the BH 
mass with less scatter than stellar velocity dispersion (Geb- 
hardt et al. 2003, but see also Marconi & Hunt 2003, who find 
a comparable scatter for the relation between BH mass and 
host bulge mass.). 

Fundamental to the understanding of the M-a relation is 
the measurement of the relation's intrinsic or cosmic scatter, 
as distinct from scatter due to measurement errors. Indeed, the 
fact that there is a relation between BH mass and stellar veloc- 
ity dispersion is not surprising, but the scatter is remarkably 
small, estimated by Tremaine et al. (2002) to be no larger than 
0.25-0.3 dex. Novak et al. (2006) carried out an extensive in- 
vestigation of the residuals from proposed M-a relations and 
variants. Their work highlights the critical role of understand- 
ing measurement errors in assessing the scatter of the various 
relations. 

The magnitude of the intrinsic scatter is extremely impor- 
tant for several reasons. First, the range of BH masses in 
galaxies of a given velocity dispersion or bulge luminosity 
constrains BH formation and evolution theories. For the past 
several years, many theories of BH formation and galaxy evo- 
lution have used the M-a relation either as a starting point for 
further work or as a prediction of the theory (e.g.. Silk & Rees 
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1998; Burkert & Silk 2001; Adams et al. 2001, 2003); for a 
review, see Richstone (2004). A further test of such theories 
is whether they can reproduce the observed cosmic scatter in 
the relation. Some predictions that are testable in principle al- 
ready exist; for example, Volonteri (2007) predicts that there 
should be an increased intrinsic scatter in low-mass galaxies 
because BHs are ejected by asymmetric gravitational wave 
emission and low-mass spheroids have lower escape veloci- 
ties. 

Understanding the scatter in the M-a relation is also essen- 
tial for estimating the space density of the most massive BHs 
in the local universe. One of the most useful aspects of the M- 
a relation is that it allows one to estimate a galaxy's central 
BH mass from the more easily measured velocity dispersion. 
Because of the steep decline in number density of galaxies 
having high velocity dispersion (Sheth et al. 2003; Bemardi 
et al. 2006; Lauer et al. 2007b), the majority of the extremely 
large BHs will reside in galaxies with moderate velocity dis- 
persions that happen to contain BHs that are overmassive for 
the given velocity dispersion (Yu & Tremaine 2002; Marconi 
et al. 2004; Lauer et al. 2007b). Knowing the magnitude of 
the intrinsic scatter is thus required to find the density of the 
most massive BHs. For example, the number density of BHs 
with M > 10^° Mq is ^ 3 Gpc~^ if the intrinsic scatter is 0. 15 
dex and ^ 30 Gpc~^ if the intrinsic scatter is 0.30 dex (Lauer 
et al. 2007b). 

Both the magnitude of the intrinsic scatter and its distri- 
bution (e.g., normal or log-normal in mass) are also impor- 
tant to know for studies of the evolution of the M-a relation 
(e.g., Treu et al. 2004, 2007; Hopkins et al. 2006b; Peng et al. 
2006; Shen et al. 2007, 2008; Vestergaard et al. 2008). Lauer 
et al. (2007a) showed that there is a bias when comparing BH 
masses derived from observations of inactive galaxies at low 
redshifts to BH masses from active galaxies at higher redshift. 
The bias arises because the sample of nearby galaxies mea- 
sures the distribution of BH masses for a given host veloc- 
ity dispersion or luminosity, whereas the sample from high- 
redshift galaxies tends to measure the distribution of the host 
luminosity or host velocity dispersion for a given BH mass. 
Lauer et al. (2007a) found that the bias in the inferred log- 
arithmic mass scales as the square of the intrinsic scatter in 
logarithmic mass. In order to account for this bias correctly, 
not only the magnitude but also the distribution of the devia- 
tions from the M-a relation is needed. 

Given the importance of the intrinsic scatter, we focus on a 
detailed examination of the scatter in the M-a and M-L rela- 
tions. This paper addresses the two most fundamental ques- 
tions regarding the scatter: (1) What is the magnitude of the 
intrinsic scatter? and (2) What is the shape of the distribution 
of the intrinsic residuals? A central part of this paper is ex- 
amination of the intrinsic scatter from galaxies lying in a nar- 
row range of velocity dispersion. This sample of galaxies in- 
cludes measurements from a companion paper (Giiltekin et al. 
2009) that presents BH mass measurements for five galaxies 
selected to fall within a narrow range of velocity dispersion 
(180 km s~^ < a < 220 km s~'). By focusing on these 
galaxies and the others in this narrow range of velocity disper- 
sion, we may study the distribution of BH masses at a given 
value of a. 

We also combine these new mass measurements with BH 
mass measurements and upper limits from the literature to 

provide new estimates of the parameters of the M-o and M- 
L relations and their intrinsic scatter. We discuss our sample 



selection in § 2. Results of fits to M-o and M-L as well as an 
analysis of their scatter are presented in § 3. We also demon- 
strate in § 4 that a bias is incurred by selecting a sample of 
BHs based on the resolution of the sphere of influence; this is 
our reason for including all reliable BH masses, regardless of 
resolution. The implications of our results are discussed in § 5 
and summarized in § 6. The details of our likelihood method 
are described in Appendix A, and tests of various error distri- 
butions and subsample selections are provided in Appendix B. 

2. SAMPLE OF BLACK HOLE MASSES 

Our analysis of the M-o and M-L relations uses a Mbh 
sample comprising the entire corpus of measurements pub- 
lished before 20 Nov 2008, augmented with the four new 
measurements (and one additional upper limit) presented 
by Giiltekin et al. (2009). We surveyed the literature to 
make a list of dynamically detected central BHs, starting 
with the compilations of Tremaine et al. (2002), Marconi & 
Hunt (2003), Ferrarese & Ford (2005), and Graham (2008). 
We include only direct dynamical measurements, and since 
reverberation-mapping measurements (e.g., Peterson et al. 
2004) are normalized to the M-o relation (Onken et al. 2004), 
they are not included here. Tables 5, 6, and 7 list all galax- 
ies with dynamical BH measurements, with masses scaled 
from the original publications to our preferred distances (also 
listed) assuming Mbh °= D and Hq = 70 km s~'Mpc~^ We 
also include upper limits to BH masses, listed in Table 6. 
Most upper limits are not restrictive, in the sense that they 
are consistent with a wide range of values for M-o intercept, 
slope and intrinsic scatter. Our fits only use BH masses from 
elhptical galaxies or from galaxies with classical bulges or 
pseudobulges (Kormendy & Kennicutt 2004). Masses from 
some previous papers by our "Nuker team" (Gebhardt et al. 
2000c, 2003) were 9% smaller than the correct values used 
here due to an error in units conversion. To correct this, the 
published values of Mbh, Miow, and Mhigh should be multi- 
plied by 1.099. This correction is small compared to the mea- 
surement uncertainties, and it is indicated in Table 5 by a su- 
perscript beside the galaxies to which it has been applied. 

It has been argued (e.g., Ferrarese & Ford 2005) that BH 
mass determinations are unreliable if the kinematic observa- 
tions do not resolve the sphere of influence of radius /?infl = 
GMbhCJ^- In particular the argument generally requires that 
the full width at half-maximum intensity (FWHM) spatial res- 
olution of the kinematic observations, i/res> satisfy /?infl/i^res > 
I to avoid bias in the BH mass determination. We do not 
agree with this argument: our tests (Gebhardt et al. 2003; Ko- 
rmendy 2004) show that the smaller values of ^infi/t^res lead 
to large error bars but not systematic bias. We expand on the 
misconception that under-resolved spheres of influence lead 
to biased BH mass determinations in § 4.1. On the other 
hand, we show in § 4 that excluding measurements on the 
basis of /?infl/<ircs leads to systematic bias, not in the individ- 
ual BH masses, but in the estimates of the parameters of the 
M-o and M-L relations. For these reasons we do not adopt 
any resolution-based cutoff in measurements included in our 
sample. We also investigate each BH mass detection individ- 
ually and decide whether any systematic uncertainties suggest 
that it should be omitted from our sample. In Table 7 we list 
published BH masses that we omit from our fitting sample as 
well as the reasons. In general, galaxies are omitted because 
(1) the authors of the study, themselves, expressed doubts of 
their models' ability to securely determine the mass, (2) there 
is no quantitative analysis of how weU their model fits the data 
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(e.g., the value of reduced for the best fit), or (3) the pro- 
vided quantitative analysis of goodness of fit is poor Note 
that one of the major conclusions of this paper is that the in- 
trinsic scatter of the M-a relation is larger than most previous 
studies have been found. If we did not omit the galaxies in 
Table 7, the scatter would increase further. After removing 
these BH masses, we are left with a sample including upper 
Umits (SU, see Table 2 for abbreviations used to denote our 
samples) and a sample excluding upper limits (S). 

In Appendix B we also consider an additional sample. Fol- 
lowing the suggestion of Ferrarese & Ford (2005) and others, 
we create a "restricted sample" (RS) of only 20 galaxies in 
which we (1) require Rinn/dres > 1-0, (2) use only detected 
BH masses (i.e., no upper limits), (3) exclude masses deemed 
suspicious by Ferrarese & Ford (2005) or by Tremaine et al. 
(2002) in their final "culled sample," (4) exclude any galaxy in 
which multiple measurements are inconsistent, and (5) make a 
subjective judgment that the quality of the mass determination 
is adequate. Because of the bias introduced by restrictions 
based on resolution, we do not reconnmend this approach, but 
we present fits to Sample RS in Appendix B. 

Membership in all samples (SU, S, RS) is given in Tables 5 
and 6. We plot the masses of the BHs as a function of velocity 
dispersion in Figure 1. 

2.1. Velocity Dispersion 

Whenever possible we use the effective velocity dispersion 
Ge as defined by 



a3 = ^ 



(1) 



where Rg is the effective radius of the galaxy and V is the ro- 
tational component of the spheroid. If this is not available, we 
use the central stellar velocity dispersion (o^) found in Hyper- 
Leda. In galaxies for which both are available, we compare 
Ge to Gc in Figure 2 and find no systematic bias to high or low 
values. Because of stellar template mismatches and possible 
errors in the determination of Rg, we impose a minimum error 
in Gg of 5%. The contribution to velocity dispersion arising 
from the relatively small rotational component of disks in spi- 
ral galaxies is unlikely to lead to large systematic errors in 
velocity dispersion. For this reason, we include spirals in our 
main sample. 

2.2. Luminosities 

We use extinction-corrected, bulge (which we obvi- 
ously interpret as total luminosity for elUpticals), V-band 
luminosities-calculated from the extinction-corrected magni- 
tudes, MO ,„,g^, using log{Lv /L,,.v) = 0.4(4.83 - M« ,^ig J. 
The choice of V band is a compromise between more widely 
available luminosities in the and V bands and less extin- 
guished ^T-band luminosities. Because bulge-disk decomposi- 
tion of spiral galaxies is fraught with difficulties and possible 
systematic errors, we do not include spirals in the Mbh-L fits. 
We do include bulge luminosities for SO galaxies for which 
we are confident in the bulge-disk decomposition because the 
disk is faint and contains little dust and few young stars. Er- 
rors in bulge-disk decomposition were estimated by examin- 
ing the range of values in the hterature. Thus the galaxies in- 
cluded in M-L fits are those which have a My ^uige value listed 
in Table 5 or 6. Most of our V-band absolute magnitudes (to- 
tal and bulge) come directly from de Vaucouleurs et al. (1991) 
as compiled by Lauer et al. (2005) and Lauer et al. (2007c). 



3. MEASUREMENT OF THE M-a AND M-L RELATIONS, AND 
THEIR INTRINSIC SCATTER 

3.1. Best Fit M-G Relation 

We are principally interested in predictors of the form 



log (M/Me ) = a + p log (o, /200 km s 



-n 



(2) 



with an intrinsic or cosmic scatter eo that is the root-mean 
square (rms) deviation in log(MBH/M0) from this relation 
(for zero measurement error). We assume for simplicity that 
£o is independent of Gg. The details of the fitting method 
are discussed in Appendix A. Based on tests in § 3.2, a log- 
normal distribution of intrinsic scatter eo is an adequate de- 
scription of the scatter. Based on tests in Appendix B, a gener- 
alized maximum-likelihood method is appropriate to use here. 
Given these, the best-fit relation for the full sample (SU) is 

(3) 

which has an intrinsic rms scatter eo = 0.44 ±0.06. These val- 
ues are very close to those obtained by Tremaine et al. (2002). 
We show in Appendix B that the choice of objects to include 
in the sample can have a significant impact on the slope, in- 
tercept, and residual scatter, but that the error distribution as- 
sumed (among those we considered) has little impact. 

The full sample may be split into subsamples: early-type 
(elliptical and SO) and late-type (spiral) galaxies, ellipticals 
and non-ellipticals, BH mass measurements by gas dynam- 
ics and other methods, low Gg and high Gg, and barred and 
non-barred galaxies. The fits are sunamarized in Table 1. In 
most cases the fits to different subsamples are consistent with 
each other. The intrinsic scatter for the ellipticals-only sam- 
ple, however, is eo = 0.31 ±0.06, which is ^2068^^ smaller 
than either the full sample (eo = 0.44 ± 0.06) or the sample 
of non-ellipticals (eo = 0.53 ± 0.10). This may reflect either 
greater unaccounted errors in BH mass and Gg measurements 
in later-type galaxies, or that ellipticals lie closer to the ridge 
Une of the M-g relation than later- type galaxies. 

3.2. Examination of the Intrinsic Scatter in M-g 

The distribution of residuals from the M-G relation is of 
practical interest as discussed in § 1. Figure 3 is a histogram 
of the residuals in log (Mbh) in the best-fit M-g relation from 
sample S. The distribution of the residuals appears consistent 
with a normal or Gaussian distribution in logarithmic mass, 
although the distribution is noisy because of the small num- 
bers. For a more direct test of normaUty we look at log(MBH) 
in galaxies with Gg between 165 and 235 km s~\ correspond- 
ing to a range in log (0^/200 kms^') from approximately 
—0.075 to 0.075. The predicted masses for the 19 galaxies 
in this narrow range differ by at most a factor of 4.3, given 
our best-fit relation. The power of having a large number of 
galaxies in a narrow range in velocity dispersion is evident 
here, as there is no need to assume a value for the slope of 
M-G or even that a power-law form is the right model. The 
only assumption required is that the ridge line of any M-G 
relation that may exist does not change substantially across 
the range of velocity dispersion. The mean of the logarithmic 
mass in solar units is 8.16, and the standard deviation is 0.45. 
The expected standard deviation in mass is 0.19, based on the 

" We use Cgg to mean 68%-confidence level so as to distinguish it from 
velocity dispersion. 
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TABLE 1 
M-a Relation for Sub samples 



Subsample 






a 


p 




^9 


Full sample 


49 


18 


8.12±0.08 


4.24 ±0.41 


0.44 ±0.06 


0.0004 ±0.018 


Early type 


38 


6 


8.22 ±0.073 


3.86±0.380 


0.35 ±0.031 


0.0145 ±0.031 


Late type 


11 


12 


7.95 ±0.286 


4.58± 1.583 


0.56±0.141 


0.0006 ±0.040 


Ellipticals 


25 


2 


8.23 ±0.084 


3.96 ±0.421 


0.31 ±0.063 


0.0006 ±0.018 


Non-ellipticals 


24 


16 


8.01 ±0.156 


4.05 ±0.831 


0.53 ±0.097 


0.0010 ±0.031 


Stars and masers 


32 


2 


8.11±0.107 


4.05 ±0.554 


0.49 ±0.075 


0.0002 ±0.021 


Gas dynamics 


17 


16 


8.16±0.122 


4.58 ±0.652 


0.35 ±0.096 


0.0036 ±0.040 


de < 200 km s"' 


25 


16 


8.07 ±0.172 


3.97 ±0.869 


0.50 ±0.091 


0.0013 ±0.031 


Oe > 200 km s~' 


24 


2 


8.12±0.158 


4.47 ±0.921 


0.35 ±0.079 


0.0026 ±0.024 


Non-barred 


41 


7 


8.19±0.087 


4.21 ±0.446 


0.43 ±0.064 


0.0006 ±0.017 


Barred 


8 


11 


7.67±0.115 


1.08± 0.751 


0.17 ±0.078 


0.1809 ±0.147 


Classical bulges 


39 


16 


8. 17 ±0.086 


4.13 ±0.434 


0.45 ±0.066 


0.0009 ±0.024 


Pseudobulges 


10 


2 


7.98±0.156 


4.49 ±0.903 


0.28 ±0.096 


0.0034 ±0.037 



Note. — Results from fits to subsamples of our full sample, based on morphological type, BH 
mass-measurement method. iVm and W,, are the number of galaxies in each group with BH mass 
measurements and upper limits, respectively. 

rms dispersion of log(Og/200 km s~^) (0.046) in this range 
times the M-o slope p; thus the variation in the ridge hne 
of the M-o relation in this sample is negligible compared to 
the intrinsic scatter. We perform an Anderson-Darling test for 
normality with unknown center and variance on this sample of 
logarithmic masses (Stephens 1974; Press et al. 1986) and find 
that the distribution is consistent with normality at better than 
the 15% level. The same sample fails an Anderson-Darling 
test for normality in Mbh — as opposed to log(MBH) — at 
the 1% level. 

The Anderson-Darling tests show that log-normal (i.e., 
Gaussian in logarithmic mass) is an acceptable description 
of the distribution of the residuals. The residuals may 
also be well represented by other distributions, which may 
be compared to a log-normal description with an odds ra- 
tio (Eq. A18). The calculated odds ratio of Gaussian to 
Lorentzian in logarithmic mass is 14.49, of Gaussian to 
double-sided exponential in logarithmic mass is 1.78, of 
Gaussian to the sum of two Gaussians in logarithmic mass 
is 1.55, and of Gaussian to a Gaussian with different standard 
deviations above and below the mean is 2.30. In these cases 
considered, the odds ratio {%ah) of the Gaussian distribution 
to the other distributions in logarithmic mass is greater than 
unity. Thus the Gaussian distribution is favored. While it 
should be noted that these tests convolve the intrinsic disper- 
sion with variance from imprecise measurements, the resid- 
uals in the M-a relation are well described by a Gaussian 
distribution in logarithmic mass. 

One component of the intrinsic scatter in both the M-a and 
M-L relations is sure to be random errors in distance. 

All BH mass measurements in our sample scale other than 
the MUky Way linearly with the distance. For these galaxies 
random errors in distance contribute directly to uncertainty 
in the mass and thus contributes to the intrinsic scatter. The 
Milky Way, which scales approximately as Z)'^ (Ghez et al. 
2008), has the uncertainty in distance incorporated in its mass 
uncertainty in our tables. The contribution of random errors 
in distance to the intrinsic scatter, however, is sure to be a 
small as the random errors in distance are typically around 
10% or 0.04 dex (e.g., Tonry et al. 2001) compared to the 



scatter of 0.44 dex. Another source of intrinsic scatter may 

be unaccounted systematic errors in Mbh measurements. In 
this paper, we assume that the measurement errors accurately 
reflect the uncertainty, but we discuss some possible cuases in 
§ 5.1. If systematics are large, then they may be an important 
contribution to the inferred intrinsic scatter. 

3.3. Log-Quadratic Fits to M-a 

We may also fit a log-quadratic function to the data as sug- 
gested by Wyithe (2006a): 



1 ( '^^ \ 

°^V200kms-i/ 



, , • (4) 
.200kms-iy'J 

For the full sample, we find a = 8.08 ± 0. 10, P = 4.47 ± 0.50, 
Y= 1.72 ± 1.71, and eo = 0.44 ±0.06. The results are con- 
sistent with log-linear (y= 0) at the lagg level, and the in- 
trinsic scatter is not significantly decreased from a log-linear 
model. The value and significance for yis similar to that found 
by Wyithe (2006b). The odds ratio (eq. A18) of a log-linear 
model to a log-quadratic model is %ib = 4.68, indicating the 
log-linear model is favored. 

3.4. Best-Fit M-L Relation 

We also look at fits to the M-L relation. Spiral galaxies, 
which present problems in determining the bulge luminos- 
ity because of the difficulty in getting a precise bulge-disk 
decomposition, are excluded. We therefore limit our fits to 
elhpticals and those SO galaxies for which we have rehable 
bulge-disk decomposition. The sample of galaxies used may 
be discerned from Tables 5 and 6 by the presence of a value in 
the column for bulge magnitude (A^}/ bulge)' which in the case 
of elhpticals is equal to the total magnitude. 

Using the same fitting method as before on all galaxies 
meeting the above criteria, we find 

logf ) = (8.95 ±0.11) + (1.11 ±0.18) log f 



(5) 

with an intrinsic scatter of En = 0.38 ± 0.09. 
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Fig. 1. — The M-a relation for galaxies with dynamical measurements. The symbol indicates the method of BH mass measurement: stellar dynamical 
(pentagrams), gas dynamical (circles), masers (asterisks). Arrows indicate SOgg upper limits to BH mass. If the SCjg limit is not available, we plot it at 3 times 
the lOfig or It 1.5 times the 2068 limits. For clarity, we only plot en'or boxes for upper limits that are close to or below the best-fit relation. The color of the error 
ellipse indicates the Hubble type of the host galaxy: elliptical (red), SO (green), and spiral (blue). The saturation of the colors in the error ellipses or boxes is 
inversely proportional to the area of the ellipse or box. Squares are galaxies that we do not include in our fit. The line is the best fit relation to the full sample: 
Mbh = 10** M ; (a/200 km s^' )'''^''. The mass uncertainty for NGC 4258 has been plotted much larger than its actual value so that it will show on this plot. 
For clarity, we omit labels of some galaxies in crowded regions. 
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(km s ') 

Fig . 2 . — Central stellar velocity dispersion 0^ from HyperLEDA as a func- 
tion of effective stellar velocity dispersion (5^ (eq. 1). There is no systematic 
bias to high or low values. The dotted line shows c^. = Oj,. The solid line 
is the best-fit regression, Oc = (-2.7± 1.4 km s"') + (1.00±0.02)af. The 
dashed line is a fit that includes uncertainties in both variables as well as an 
intrinsic scatter: ac = (-0.9±7.2 km s"' ) + (0.95 ± 0.09)a(, with intrinsic 
scatter of 22 ± 5 km s^' . Both fits are consistent with c^. = (5^ within their 
uncertainties. 
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Fig . 3 . — Histogram of residuals from the best-fit M-o relation in sample S. 

3.5. Examination of the Intrinsic Scatter in M-L 

The scatter determined here for the M-L relation is notably 
smaller than other studies have found, and it is consistent 
with the intrinsic scatter in M-o in early-type galaxies. Other 
studies have found similar scatter between the M-L and M- 
O relations (e.g., Marconi & Hunt 2003). Figure 5 shows a 
histogram of the residuals from the fit. We test the distribu- 
tion of 12 masses in a narrow range in luminosity (10.2 < 
log(Lv/L0 v) < 10.7) for normality and log-normality using 
an Anderson-Darling test with unknown mean and variance. 



The mean logarithmic BH mass in that range is 8.21 with stan- 
dard deviation 0.36. The expected standard deviation in loga- 
rithmic mass is 0. 16, if all the galaxies lie exactly on the ridge 
line of the M-L relation, based on standard deviation in log- 
arithmic luminosity of 0. 14. Thus the variation in the ridge 
line of the M-L relation in this sample is negligible compared 
to the intrinsic scatter The distribution of masses is consis- 
tent with a log-normal distribution, but is inconsistent with a 
normal distribution at the 1% level. 

Again, we consider several possible functional forms of the 
intrinsic scatter in logarithmic mass at constant luminosity 
and calculate odds ratios of Gaussian to Lorentzian (3.37), to 
double-sided exponential (1.10), to double Gaussian (1.36), 
and to a Gaussian with different standard deviations above 
and below the mean (1.99). Thus, normal or double-sided ex- 
ponential in logarithmic mass are equally acceptable descrip- 
tions of the intrinsic scatter. There are, however, relatively few 
mass measurements in this range, so this it is not as strong a 
test as it is for M-o. 

4. BIASES INTRODUCED WHEN CULLING THE SAMPLE BY 
BLACK HOLE SPHERE-OF-INFLUENCE RESOLUTION 

Our analysis has used all Mbh measurements in the sample 
without regard to how well the sphere of influence is resolved. 
In essence, we trust the uncertainties returned by the dynam- 
ical models to reflect the quality of the Mbh determination — 
given the assumptions — for whatever resolution was realized. 
This approach conflicts with the caveats advanced by, e.g., 
Ferrarese (2002), who concluded that Rmfi / d^es > 1 was a nec- 
essary (but not sufficient) condition to make a precise Mbh 
determination. Such concerns have motivated a number of 
groups studying the BH relations (e.g., Merritt & Ferrarese 
2001b; Ferrarese & Ford 2005) to cull the more poorly re- 
solved BHs from their samples. We show in this section that 
far from improving the accuracy of the estimated relations, 
this procedure will actually produce biased estimates of the 
intercept, slope, and intrinsic scatter of the relations. 

4.1. What Resolution of the Sphere of Influence Does and 
Does Not Mean 

We advocate use of the sphere-of-influence scale only as a 
rough guide to the needed spatial resolution of the observa- 
tions. In the arguments that follow, we will show that Mbh 
determinations from datasets with Rj^n / '^res < 1 are unbiased 
when using three-integral models. This, however, contrasts 
with claims commonly made in other works: that strict resolu- 
tion of the sphere of influence is required for credible Mbh de- 
terminations and more importantly, that Mbh determinations 
made from observations that do not resolve the sphere of in- 
fluence will be biased. Given the strong prevalence of this 
viewpoint, and prompted by comments from the referee, we 
review its development and application in the literature. We 
find, in fact, that there is little or no support for the conclusion 
that Mbh determination become increasingly biased with de- 
creasing resolution. It appears that the common but uncritical 
application of sphere-of-influence-resolution as a way to cull 
Mbh determinations cannot be justified by careful reading of 
the very works often cited in its support. 

In their review article, Ferrarese & Ford (2005) write, "All 
studies which have addressed the issue [of BH mass determi- 
nation and resolution level]. . . have concluded that resolving 
the sphere of influence is an important (although not suffi- 
cient) factor: not resolving [Ri„f[] can lead to systematic errors 
on Mbh or even spurious detections," and cite the following: 
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Fig. 4. — The M-L relation for galaxies with dynamical measurements. The symbol indicates the method of BH mass measurement: stellar dynamical 
(pentagrams) and gas dynamical (circles). Arrows indicate upper limits for BH mass. Squares are galaxies that we omitted from the fit. The color of the error 
ellipse indicates the Hubble type of the host galaxy [elliptical (red) and SO (green)] and the saturation of the color is inversely proportional to the area of the 
ellipse. The line is the best-fit relation for the sample without upper limits: Mbh = 10*" M,t (Lv /lO" L,t, y)' " . 



Ferrarese & Memtt (2000); Merritt & Ferrarese (2001b,a); 
Graham et al. (2001); Ferrarese (2002); Marconi & Hunt 
(2003). We consider each of these in turn. 

Ferrarese & Merritt (2000) found that the ground-based 
Mbh measurements by Magorrian et al. (1998) were higher 
for fixed velocity dispersion than the predictions of their em- 
pirical M-a relation and judged them to be therefore biased. 
The discrepancy with their M-a relation increased with in- 
creasing distance. While discrepancy with the M-a relation 
is not a justifiable reason for excluding Mbh measurements 
from the relation (the argument is circular), the masses from 
Magorrian et al. (1998) were, in fact, biased to high values 



by roughly a factor of 3. The reason for the bias, however, 
is that they came from two-integral, isotropic, axisymmetric 
models, not because they were more poorly resolved (Merritt 
& Ferrarese 2001a; Gebhai'dt et al. 2003). 

Merritt & Ferrarese (2001b) present similar arguments as 
do Merritt & Ferrarese (2001a) who also go on to describe 
the reason two-integral models yield masses that are biased 
somewhat high. Merritt & Ferrarese (2001a) do mention that 
three-integral models provide an increased space of orbits that 
can lead to an increased range in acceptable black hole masses 
but not that they are biased, resolved or not. 

Graham et al. (2001) find decreased scatter in the relation 
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Fig. 5. — Histogram of residuals from best-fit M-L relation. 



between Mbh and concentration index when removing more 
poorly resolved galaxies in addition to another galaxy. Again, 
using a provisional relation to exclude potential data points 
from the relation is not the same as the identification of a bias. 

Ferrarese (2002) does not present new arguments or studies 
but repeats the arguments just described. 

Marconi & Hunt (2003) say nothing about a bias, only 
that they used /?infi/<ires as a criterion for determining reliable 
masses and that they find a smaller scatter in the /T-band M-L 
relation when excluding lower resolution observations. 

Valluri et al. (2004) are frequently cited as having shown 
that resolving the sphere of influence is necessary for accu- 
rate Mbh determination, but, in fact, they do not make such a 
claim. The most important result from their work was to show 
that when using too few orbits, spurious black hole mass de- 
terminations may appear. When enough orbits are used, Val- 
luri et al. (2004) found that a wide range of values for Mbh 
are acceptable with their synthetic data set, but this turns out 
to be caused by the lack of real measurement noise in the data 
set (Magorrian 2006). Even so, the input Mbh in their sim- 
ulations was always within their range of acceptable values, 
i.e., not biased. The bias that Valluri et al. (2004) do identify 
is what results from over-regularizing (i.e., requiring smooth- 
ness in orbit solutions) not from under-resolving the sphere of 
influence. 

The degree of resolution of the sphere of influence of the 
BH cannot tell the complete story of the reliability of the BH 
mass determination. In the quantity Rinn/dyes there is encoded 
no information about, e.g., the spectroscopic resolution of 
the data. To illustrate this, we perform a simple experiment. 
We take surface brightness profiles, /(r), from two galaxies 
with velocity dispersion Oe ~ 200 km s^\ one a core galaxy 
(NGC 3607) and one a power-law galaxy (NGC 4026), and 
deproject to get luminosity densities. We assume a constant 
mass-to-light ratio of T = 4 and calculate the projected veloc- 
ity dispersion profile, Oo(r) from a spherical, isotropic model, 
assuming a seeing of 0"1. We repeat this process but this 
time assuming there is a black hole with mass Mbh = lO^M v, 
at the center to get aBH('")- We then calculate the difference 
in projected velocity dispersions from the two models. The 
ability to discriminate between these two models depends on 
the measurement error of velocity dispersion. So. We calcu- 



late the significance of the difference between these models 
[Obh('") — Oo('")]/§o for three values of the measurement er- 
ror 8o = 5, 10, and 20 km s^\ which roughly correspond to 
signal-to-noise ratios from HST STIS data of more than 50, 
40 to 50, and 35, respectively, depending on the details of 
the observation. We plot these curves in Figure 6 as solid 
lines. Two conclusions are immediately obvious from these 
curves: (1) It is possible to discern a significant difference 
between the two profiles outside of r/R-mn = 2 (correspond- 
ing to /?infl /fires < 1), and (2) the significance of the differ- 
ence depends on the measurement errors. An actual observa- 
tion would amplify the difference because, rather than com- 
paring the two velocity dispersion profiles at one location in 
radius, it would integrate the velocity dispersion within a reso- 
lution bin. To illustrate this we calculate Pbh('') — J^o('')]/8o, 
where 



IoKr)rdr 



(6) 



and similarly for Eq. We plot these as dashed lines in Fig- 
ure 6 for each of the representative measurement errors. These 
curves show that it is possible to determine the presence of a 
BH when your data do not resolve the sphere of influence, 
depending on the quality of the data. 

It is further evident from the figure that for a given measure- 
ment error in velocity dispersion, the ability to discriminate 
between different BH masses decreases further away from 
the center, i.e., errors in BH mass increase as /?infl/iires de- 
creases. This fact has been seen several times before in the 
literature. Gebhardt et al. (2003) modeled 12 galaxies using 
(1) ground-based data only and (2) ground-based data com- 
bined with space-based data. In every case, the models re- 
cover larger uncertainties for the data sets with only ground- 
based data. A separate study, by Kormendy (2004), consid- 
ered the many mass measurements of the BH in M32 using 
a variety of modeling techniques and data that range in level 
of resolution by over an order of magnitude. The study re- 
vealed that the error bars in the Mbh determination increased 
as Rinn/dies decreased but that the Mbh values were not bi- 
ased at any resolution. In fact, the value for Mbh found from 
the highest resolution data was consistent with all of the lower 
resolution results. This held true despite the fact that the val- 
ues considered were derived from widely varying modeling 
methods. 

When looking at the data in Table 5, it is tempting to look 
for a correlation between the size of the error in BH mass and 
Rinn/dres, but the size of the error depends on many factors. 
Among the factors are differences in method of BH mass de- 
termination, differences in codes within similar methods, and 
varying quality data among the different galaxies. When con- 
sidering only stellar dynamical mass determinations by the 
code described by Siopis et al. (2008), there is a weak yet sig- 
nificant trend, as seen in Figure 7, which plots Mhigh/Miow as 
a function of Riaa/dies- This trend is weak because of other 
factors that determine the precision of the measurement. 

We stress again that larger error bars do not imply a bias 
in the value of BH mass recovered from observations with 
smaller Rmn/dres- For example, Gebhardt et al. (2003) show 
that there is no significant bias in their figure 9, which plots 
BH masses obtaimed from ground-based data alone as a func- 
tion of masses from the same ground-based data combined 
with space-based data. The line of equality in their figure goes 
through 10 of the 12 error bars. In 4 of the 12 galaxies, the 
difference between the best-fit values is less than 25%, and in 
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Fig. 6. — Plots of significance of tlie difference between velocity dispersions of a galaxy with a BH and one without. The solid lines show (cbh — Oo)/8a as a 
function of radial distance in units of Rinfl for three different values of measurement error in velocity dispersion. The dashed curves show the difference between 
light-weighted integrated velocity dispersion profiles. The left panel is derived from NGC 3607, a galaxy with a core surface brightness profile, and the right 
panel is derived from NGC 4026, a galaxy with a power-law surface brightness profile. It is clear that (1) it is possible to discern the presence of a BH outside of 
the sphere of influence and (2) the ability to discern the presence of a BH depends on the measurement error. 

synthetic M-a data set. Each synthetic data set consists of a 
sample of 40 galaxies, uniformly distributed in volume out to 
a distance of 30 Mpc. Each galaxy is given a velocity disper- 
sion from a normal distribution in log(o/200 km s^^) cen- 
tered at with standard deviation 0.2. Each galaxy is given 
a BH mass from an M-a relation with a = 8, P = 4.0, and 
log-normal intrinsic scatter with eo = 0.3 dex. The BH's log- 
arithmic mass is measured with a normally distributed mea- 
surement error of 0.2 dex and the velocity dispersion has a 5% 
error. Since each galaxy has a distance, a BH mass, and a ve- 
locity dispersion, we calculate Ri^n = GMbhO^ and assume 
the galaxy to be observed with an instrument with resolution 
of c/i-es = 0"1. This process is repeated for 10^ realizations for 
each of the choices in assembling the sample described below. 

We first tried fitting for the parameters of the M-a relation 
using these simulated data sets and the same fitting proce- 
dures that we applied earlier to the actual data. We success- 
fully recovered the input parameters with only a slight bias 
to low intrinsic scatter (recovering e — 0.27 instead of the in- 
put value of 0.3 dex). The samples were then culled in two 
ways: (1) we removed galaxies that fell below a given cut- 
off in Rina/dres, where Ri„f[ is computed from the measured 
BH mass and velocity dispersion; and (2) we removed galax- 
ies that fell below a given cutoff in Rmn/dres, where Rian is 
computed from the BH mass predicted by the M-a relation. 
Since the Galaxy plays an important role in the observed sam- 
ple (because of its position near the low end of the range of 
observed o^. and the small uncertainty in the measured mass), 
we also augmented some of the samples by adding one galaxy 
with the actual measured values of the Galaxy's BH mass and 
velocity dispersion with their corresponding uncertainties. 

4.3. Culling the Simulated Sample Based on Traditional 
Sphere of Influence 

First, we use various cuts in Rm^/dres to eliminate galaxies 
from the simulated sample. The results from this simulation 
are shown in Figure 8. Without the Galaxy, the trends are sim- 
ple and monotonic: increasing the /?infl/<ires cutoff increases 
the value of the intercept and decreases the values of the slope 
and scatter. Because our sample RS (see Appendix B) omits 
the Galaxy, comparison of the orange curve to the black curve 
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Fig. 7. — Plot of size of errors in BH mass, calculated as Mi,igh/Mio„, as a 
function of i?i„fl /rfrcs for stellar dynamical measurements using the code de- 
scribed in Siopis et al. (2008). There is a weak but significant trend indicating 
that lower resolution contributes to larger errors in BH mass, but it is not the 
only factor in determining the uncertainty. 

9 of the 12 it is less than 50%. While 9 of the 12 have smaller 
values of Mbh for ground-only datasets (actually the oppo- 
site bias to that seen in the Magorrian et al. (1998) models, 
which is additional evidence that the limited resolution does 
not produce a systematic bias), for 12 independent trials of 
an event with 50% chance of success, there will be 9 or more 
successes 7.3% of the time. This is entirely consistent with 
no bias. So while increased spatial resolution will always im- 
prove the precision of BH measurements, this trend appears 
to be correctly reflected in the larger error bars that emerge 
from the modeling procedure. 

4.2. Simulated Sample of Galaxies 

We demonstrate the effects of such a culled sample selec- 
tion with a series of simple Monte Carlo experiments on a 
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in Figure 8 approximates comparison of the fit values in our 
RS sample to our S sample. They show the same trends in 
intercept, slope, and scatter. 

The reasons for these biases can be easily seen in Figure 9a, 
which illustrates the effects by plotting a sample of 500 galax- 
ies randomly generated in the same way and plotting dif- 
ferent ranges in Rinn/dve^ with different symbols. Because 
/Jinfl w MbhCT^^ and Mbh ~ ci^, lines in constant Rmn/dies 
tend to fall on lines of Mbh ~ cr^~'^. (This argument ne- 
glects the different distances of the galaxies, but since the 
Monte Carlo simulation is uniformly distributed in volume, 
most galaxies are near the outer edge of the volume.) Since 
the synthesized data set uses P = 4.0, each subsample has a 
decreased scatter because the cuts fall along lines with slope 
P — 2 = 2.0 and take out the bottom portion of the scatter. The 
slopes are similarly biased to low values because the points 
have been removed from systematically smaller masses and 
systematically smaller velocity dispersions. The intercept in- 
creases as each subsample's mean increases. 

When the Galaxy is included (Fig. 8b), the magnitude of the 
biases decreases. The intercepts are still biased to high values. 
The slopes are biased to low values except those from the most 
restrictive cutoff, which are biased to a shghtly higher value. 
The scatters are biased to a smaller value; a significant number 
are consistent with eo = 0. The bias when using older, lower 
values for our Galaxy's BH mass (Ghez et al. 2005) would 
have produced more extreme changes in slopes. 

4.4. Culling the Simulated Sample Based on Sphere of 

Influence from Expected Black Hole Mass 

For the second culling strategy, we adopt a different def- 
inition of /?infl that does not lead to the extreme biases seen 
above. Because the traditional definition of /?infl depends on 
the measured BH mass, it is sensitive to whether the BH 
is overmassive or undermassive for a given velocity disper- 
sion. Thus, using /Jing as a sample criterion necessarily bi- 
ases any measurement in scatter. Instead, we offer an alterna- 
tive definition based on the expected value of Mbh from the 
M-a relation. That is, we replace culUng based on /?infl/i/res 
by culling based on /?exp/<ires. where /?exp = GMpredCJ"^ and 
Mpred = 10"(a/200 km s~')P M®. Using this as a selection 
criterion for finding the M-a relation requires only a simple 
and rapidly convergent iterative procedure. It is possible that 
this procedure will not always converge to a unique answer, 
but this is unlikely to happen if the initial guess is chosen 
close to the actual relation. Adopting this strategy substan- 
tially decreases the bias in the fitting parameters (Fig. 10) but 
dramatically increases their uncertainties since the range of a 
used in the fit is reduced. The reason for this is seen in the 
right panel of Figure 9. Instead of falling along lines of slope 
2.0, the cuts tend to fall on Unes of constant a. 

4.5. Culling Based on Traditional Sphere of Influence Using 
Observed Sample with Simulated Masses 

The above synthetic datasets represent an ideaUzed scenario 
in which the galaxies come from a volume-limited sample. 
The real sample of BH mass measurements, however, is not 
uniformly distributed in volume. As a final experiment on 
a synthetic data set, we use the observed sample of velocity 
dispersions, errors in velocity dispersion, distances, errors in 
mass measurements, and instrumental resolution. The masses 
are synthesized from an M-a model as before and censored 
according to /?mfl/<ires as calculated from the synthesized BH 



mass. There is no need to artificially include the Galaxy in 
these simulations since it is used along with all of the other 
data. The resulting distributions of fit parameters are pre- 
sented in Figure 11. Again, using the entire sample recovers 
the input parameters. Cuts in Rian/dtes show the same bias to 
high values of intercept as well as reducing the efficiency of 
the estimates. The slope is biased to high values. The intrin- 
sic scatter is similarly biased to low values, with a noticeable 
number of realizations recovering zero intrinsic scatter. We 
plot 5 realizations and indicate their level of resolution in Fig- 
ure 12. The same trends as seen in the idealized simulations 
are present. The trends in parameter estimation from culled 
samples seen here may explain why Ferrarese & Ford (2005) 
find a higher intercept (a = 8.22 ±0.06), a much higher slope 
(P = 4.86 ± 0.43), and a "negligible" scatter with their sample 
restricted to7?infi/t/res > 1- 

The systematic effects of /?infl/<ircs cutoffs can also be seen 
in Figure 13, which plots the residuals to our M-a fit as a 
function of a^. The lower resolution objects tend to be found 
toward the bottom and right of the plot, and the most highly 
resolved galaxies clearly show a residual trend in a^ even 
though the sample as a whole shows no such trend. Hence, 
eliminating reliable BH mass measurements (in the sense that 
the error bars accurately reflect the uncertainty in the mea- 
surement, even if they are large) by selecting on /?infl/c?res will 
bias the resulting fits. This bias could be mitigated by us- 
ing /?exp rather than 7?infl in the culling criterion but then the 
uncertainties in the fit become much larger. Thus, our rec- 
ommendation is to include BH masses without regard to the 

value of /?iiifl /fires- 

5. DISCUSSION 

5.1. Systematic Uncertainties in Measurements of BH Mass 

Our measurements of the intrinsic scatter in the M-a and 
M-L relations assume that the stated errors in Mbh accurately 
reflect the measurement uncertainties. If unaccounted system- 
atic errors are large, however, they will significantly increase 
the inferred intrinsic scatter. We list some potential sources of 
systematic errors in Mbh measurements. 

• Models in which the value of i smuch larger than 
the number of parameters are inadequate. We have re- 
moved Mbh measurements from our sample where this 
is a problem. 

• For more massive galaxies, the derived stellar mass-to- 
light ratio (Th,) may be wrong due to an un-modeled 
contribution of the dark matter halo to the stellar kine- 
matics at the large radii used to constrain the stellar 
population. If the dark matter halo is less centrally 
concentrated than the stellar mass distribution, as is ex- 
pected from cosmological simulations of adiabatic col- 
lapse (Gnedin et al. 2004), then the intrinsic would 
be lower than inferred, implying larger than estimated 
Mbh values (Gebhardt & Thomas 2009). 

• If the BH mass is measured by gas kinematics, uncer- 
tainty in inclination and non-gravitational forces may 
cause systematic errors (Ho et al. 2002). 

• Galaxies that have recently merged such as Cen A may 
not be in the equilibrium state that is assumed by stellar 
kinematic models. 
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Fig. 8. — Histograms of fit parameters from Monte Carlo simulations of 10^ realizations of 40 galaxies with BH masses derived from a model with a = 8, 
p = 4.0, and Eq = 0.3 with rms measurement errors in log mass of 0.2. The right-hand panels differ from the left panels in that one of the 40 galaxies generated 
by the Monte Carlo procedure is replaced by a galaxy having the mass and velocity dispersion of the Milky Way. The black curves show results when all galaxies 
are included in the fits. These curves show that our fit method is able to recover the original parameters with a slight bias to lower intrinsic scatter. The remaining 
curves come from imposing a cutoff in Rinfl/'^res and fitting the remaining galaxies from the original 40. The colors correspond to a cutoff value of 1.0 (red), 2.0 
(orange), 3.0 (green), and 4.0 (blue). The biases in the left panels can be summarized as follows: (1) the intercept a, increases monotonically with increasing 
cutoff in i?infl /^i'lcs ; (2) the slope p decreases monotonically with increasing cutoff; and (3) the intrinsic scatter Eq decreases monotonically with increasing cutoff, 
with a large number of samples consistent with Eq = 0. The orange curve (Sinti/'^res > 1) may be compared with our fits to sample RS, which does not contain the 
Galaxy. If our Galaxy is included (right panels), the biases are (1) the intercept increases monotonically with increasing cutoff in Sjnfl/t/res, though the presence 
of the Galaxy mitigates this somewhat; (2) the slope is biased to slightly lower values for the first three cutoffs and to a slightly higher value for the highest cutoff; 
and (3) the inferred intrinsic scatter is biased to a lower value, with a significant number of runs consistent with Eq = 0. 



• There are galaxies with AGN, whose central kinematics 
are hard to measure. This is systematic noise since it 
depends on how one handles the AGN subtraction. 

• Uncertain deprojection of the gas density profiles may 
be a source of systematics. 

5.2. The Value of the Slope in the M-O Relation 

The slope we report here (p — 4.24 ±0.41) is noticeably, 
though not significantly, steeper than the slope obtained by 
Tremaine et al. (2002) (p = 4.00 ±0.31). For comparable as- 
sumptions, however, this difference disappears. The most di- 
rect comparison is to apply the symmetric least-squares fitting 
method of Tremaine et al. (2002, see Table 3) to our sample 
without upper limits, which yields P = 4.02 ± 0.37, compared 
to p ^ 4.00 ±0.31 found by Tremaine et al. (2002). The two 
slopes are consistent with each other. 

5.3. The Value of the Scatter in the M-O Relation 



The value of the intrinsic scatter obtained here (eo = 
0.44 ± 0.06) is significantly larger than the value estimated 
by Tremaine et al. (2002), which was no more than Eq = 0.25 
to 0.3. Since the primary interest of the M-a relation derives 
from its tightness, a conclusion that the scatter is larger than 
previously believed warrants attention. Tremaine et al. (2002) 
used both a different fitting method and a different sample. 
We examine both of these, finding that the sample differences, 
rather than the fitting method or updated Mbh measurements 
for a given galaxy, are the cause of the different scatter mea- 
surements. 

We fitted several different samples: (1) the original 
Tremaine et al. (2002) sample, (2) the Tremaine et al. (2002) 
sample updated to the values used in this paper, (3) sample S 
from this paper, and (4) galaxies in sample S that are not in 
Tremaine et al. (2002). Fits were done using both the maxi- 
mum likelihood method developed for this paper and the sym- 
metric least-squares method used by Tremaine et al. (2002). 
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Fig. 9. — Monte Carlo sample of 500 galaxies with BH mass generated from an M-a relation with a = 8 and p = 4.0 and a log-normal scatter of 0.3 dex with 
measurement eiTors of 0.2 dex. The ridge line of the M-0 relation is drawn as a black line. The Galaxy is plotted as a pentagram. Both panels use different 
colors and symbols for different resolution levels. The left panel defines Rjnfl = GMbh/O" with Mbh equal to the measured BH mass; and the right panel 
uses our alternate definition using the expected mass of the BH based on the M-a relation. For each definition, the symbols are Si„(]/(/fi;s < 10 (red squares), 
1.0 < Sjnfl/t/ies < 2.0 (orange diamonds), 2.0 < Rinfl/'^res < 4.0 (green triangles), and Rinfl/'^res > 4.0 (blue circles). In the left panel, fitting the combined 
subsamples with Rinfl/^^rcs exceeding a given value yields a shallower slope and a smaller intercept than the underlying M-a relation.. The reason for the bias 
in slope is that cuts in Rjnfl tend to fall along lines of Mbh a^^^ (since i?i„fl « MbhO^ and Mbh o^)- This is illustrated by the dashed line of slope 2.0. 
Cuts using the expected R\^f[ (right panel) do not suffer from this bias because the cuts in R\^f[ / d^^t, with this definition tend to run at constant a. The bias they 
introduce is far less than with the traditional definition. The uncertainties in the parameters when fitting, however, are much larger due to the decreased dynamic 
range in a. 



We use our sample without upper limits to make the most di- 
rect comparison between the methods. The intrinsic scatter 
that we find in sample S is almost identical to that in sam- 
ple SU. When used on the same sample, there was never a 
significant difference between the two methods, which leads 
us to quote only the results from the maximum likelihood 
method of this paper The original Tremaine et al. (2002) 
sample yields an intrinsic scatter of e = 0.31 ±0.06. Up- 
dating the values in the same sample increases the intrinsic 
scatter to e = 0.35 ± 0.06. The intrinsic scatter in sample S 
is £ = 0.43 ± 0.06, and the intrinsic scatter from just the 19 
galaxies in sample S not found in Tremaine et al. (2002) is 
e = 0.51 ±0.11. Since the two methods used do give consis- 
tent results, they cannot be the source of the difference. 

Next we examine the difference between the original 
Tremaine et al. (2002) sample and the new galaxies. 
Kolmogorov-Smirnov tests on the distributions of mass, ve- 
locity dispersion, and error in logarithmic mass cannot rule 
out the null hypothesis that the two samples came from the 
same parent distribution at better than 50% confidence. 

A chi-square test on the distribution of morphological 
types, however, reveals that the probability that the two sam- 
ples in galaxy type come from the same parent distribution is 
only approximately 20%. This is far from conclusive, but it 
is noteworthy that the more recent Mbh measurements con- 
sist of a relatively larger number of spirals. In § 3.1, we find 
that the intrinsic scatter in the population of late-type galaxies 
may be larger than the scatter in early-type galaxies. This may 
be a result of unrecognized systematic effects in Mbh mea- 
surement of spirals or because the scatter is actually larger. 
In contrast, we find consistent intrinsic scatter estimates, at 
about 0.3, when fitting only early-type galaxies from the two 
samples. 



In this context, Hu (2008) has noted that pseudobulges ap- 
pear to host relatively smaller black holes than do "classical" 
bulges. We have included both bulge types in the analysis, and 
the differences in the relative contributions of the two types 
between various samples may enhance the intrinsic scatter. 

One of the spiral galaxies that was not included in the 
Tremaine et al. (2002) samples is Circinus, the largest out- 
lier in our sample. When we exclude Circinus from sample S, 
the scatter reduces to e = 0.36 ± 0.05, in near agreement with 
Tremaine et al. (2002). While it is disturbing that a single 
galaxy can cause such a large change in our scatter estimates, 
Circinus is an extreme outlier, a factor of ^ 30 below the M- 
O ridge line. In order for any M-a model to explain such an 
outlier that does not have large measurement errors, a large 
intrinsic scatter is needed. 

We summarize the reasons for the difference in intrinsic 
scatter measurements as follows. 

• The difference in fitting methodology is not the source 
of the difference in intrinsic scatter estimates. 

• The difference in samples is the source of the differ- 
ence. 

• Updating the Mbh measurements of galaxies in 
Tremaine et al. (2002) slightly increases the intrinsic 
scatter estimate. 

• The galaxies in the Tremaine et al. (2002) sample and 
the new galaxies in this sample possibly have different 
distributions of Hubble types. If they do, the increased 
fraction of spirals in the newly added galaxies may be 
a source of the increased scatter. The underlying cause 
may be either because (1) spiral galaxies are suscep- 
tible to more systematic error or (2) there is a larger 
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Fig. 10. — Histograms of fit parameters from Monte Carlo simulations as in Figure 8, except the cuts are based on i?exp/^'ies as defined in § 4.4. The colors 
correspond to a cutoff value of 1.0 (red), 2.0 (orange), 3.0 (green), and 4.0 (blue). The biases in the left panels are minimal, but the efficiency of recovering the 
original parameters is greatly reduced since the dynamic range in a is greatly reduced. In the right panels, the biases incuired from these cutoffs are small, and 
the efficiency of recovering the original parameters is greatly improved compared to the left panels since the dynamic range in a always includes the Galaxy. 
There is, however, still a tendency to infer an intrinsic scatter consistent with Eq = when making a large cut. 



intrinsic scatter in Mbh in spiral galaxies, as might oc- 
cur if there is a significant difference in the population 
of black holes hosted by pseudobulge versus classical 
bulges. 

• When considering only early- type galaxies, the scatter 
values of the samples are consistent with each other and 
are close to the value quoted by Tremaine et al. (2002). 

• The inclusion of Circinus makes a large difference in 
the scatter estimates. 

5.4. The Role of Upper Limits 

Upper limits provide information on the population of BH 
masses. Weak upper limits provide little information, but 
when the limit is comparable to detected BH masses at a sim- 
ilar velocity dispersion or luminosity, it can provide a valu- 
able constraint on the fit. Thus, excluding upper limits from 
the fit excludes potentially useful information. Since restric- 
tive upper limits will lie at lower masses, excluding them may 
also introduce a bias, especially in estimates of the intrinsic 
scatter and intercept. For these reasons, we view our fit in- 
cluding upper limits as the most complete answer. To avoid 



being biased by upper limits measured in galaxies with no 
BH, which do not belong in either the M-a or M-L relations, 
we include the probability of a galaxy's having no BH (fg) in 
our fit. As the simplest possible form, we assume Pq is con- 
stant for all galaxies. Given that smaller galaxies may be less 
likely to form central BHs and that their smaller escape veloc- 
ities mean ejection of their BHs is more likely, this assump- 
tion may not hold. Tests using a simple linear dependence 
on velocity dispersion, — ^ + r|log (o/200 km s^^), found 
r| = 0.001 ±0.047. So assuming a constant form is unlikely 
to introduce a strong bias, and in any event the data are not 
strong enough to support a more elaborate model. Excluding 
upper limits reduces the slope, although not by a statistically 
significant amount. 

5.5. Implications for Space Density of Black Holes and 
Studies ofM-O Evolution 

The intrinsic scatter in BH mass scaling relations affects es- 
timates of the space density of BHs based on velocity disper- 
sion or luminosity functions (Yu & Tremaine 2002; Marconi 
et al. 2004; Lauer et al. 2007b). This is especially true for 
the largest BHs. The fundamental reason is that most of the 
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Fig . 12. — Same as Figure 9, except that the velocity dispersions, distances, 
Slope and errors in mass measurements are all deiived from the actual observed 

sample of galaxies. The masses are synthesized from anM-o model. 
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Scatter 

Fig. 1 1 . — Histograms of fit parameters from Monte Carlo simulations as 
in Figure 8, except that the velocity dispersions, distances, and errors in mass 
measurements are all derived from the actual observed sample of galaxies. 
The masses are synthesized from an M-a model. 



largest BHs will come from lower dispersion or lower lumi- 
nosity hosts with overmassive BHs. We reproduce the cumu- 
lative density functions based on our results in Figure 14. In 
order to calculate the density of BHs, a predictor for Mbh 
is needed, which requires a regression on Oe or Ly. Our 
maximum-likelihood method may be used to do a regression 
by simply setting the uncertainties in Og and Ly to zero. The 
resulting fits differ very little from those presented above. 

Figure 14 shows that there is a marked difference in the 
space density of BHs when assuming no scatter or our best- 
fit scatter. Regressions for both the entire sample and for 
ellipticals-only are shown in the left panel of Figure 14. Be- 
cause the largest BHs are found in elliptical galaxies, the 
ellipticals-only curve with scatter best reflects the velocity- 
dispersion-based calculation in this paper. Lauer et al. 
(2007b) argued that the more fundamental scaling relation for 
the largest BHs may be M-L. If this is the case, then the curve 
for early-type galaxies with scatter in the right-hand panel of 
Figure 14 is the appropriate function to use. The difference 
between the two curves at the high-mass end is considerable, 
and it is a reflection of the inconsistency of the M-a and M-L 
relations in this regime. Studies of the evolution of M-a are 
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Fig. 13. — Residuals in mass from our fit to M-o for the full sample. 
The colors and shapes of the subsamples correspond to how well resolved the 
BH's sphere of influence is (from low to high Rina/ dies'- red squares, orange 
diamonds, green tiiangles, and blue circles). Removing the more poorly re- 
solved objects from the sample leaves a residual trend, even though there is 
no such trend for the sample as a whole. 
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also biased due to the high scatter. 

6. CONCLUSIONS 

In this work we have presented the latest results relating 
the mass of a galaxy's central BH with the stellar velocity 
dispersion and bulge luminosity of the galaxy. 

(1) Most current M-G fit. We compiled 49 BH mass mea- 
surements and fit a new M-a relation using a maximum- 
likelihood method that includes an intrinsic scatter compo- 
nent. We also include 18 upper limits to BH masses. The best 
fit for the relation is 

(7) 

with an intrinsic rms Gaussian scatter of Eq = 0.44 ± 0.06. 

(2) Characterization of the intrinsic scatter in M-a. Study 
of the distribution of BH masses in galaxies in a narrow range 
in velocity dispersion reveals that the intrinsic scatter in the 
M-a relation is consistent with log-normal in BH mass and 
inconsistent with normal in BH mass. Further, we find no ev- 
idence that adopting a log-quadratic, rather than a log-Unear, 
relation between BH mass and velocity dispersion reduces the 
scatter. 

(3) Smaller scatter in the population of ellipticals. When 
we limit our sample to include elliptical galaxies, we find that 
the scatter decreases to eo = 0.31 ±0.06. This is ~ 2a68 
smaller than the intrinsic scatter in non-elliptical galaxies, 
Eo = 0.53 ± 0.10. If real, this difference could either be be- 
cause of larger systematic errors in BH mass determinations 
of late-type galaxies or because elliptical galaxies he closer to 
the M-a relation ridge line. 

(4) Most current M-L fit. We also fit a relation between the 
mass of the central BH and the V-band bulge luminosity of 
the host. The best-fit for the relation is 

(8) 



with an intrinsic scatter of eo = 0.38 ± 0.09. 

(5) Characterization of the scatter in M-L fit. We also find 
through study of the distribution of BH masses in galaxies in a 
narrow range of bulge luminosity that the scatter is adequately 
described by a log-normal intrinsic scatter in BH mass and is 
inconsistent with a distribution normal in mass. 

(6) Identification of bias in samples that use Rmfl/dres a 
selection criterion. Selecting a sample for fitting an M-a rela- 
tion with a minimum value of Rmn/drts causes a bias that leads 
to an overestimate of the intercept, an overestimate estimate 
of the slope, and an incorrect estimate of the intrinsic scat- 
ter. The bias arises because Rmn ~ GMbh^^^^ and Mbh ~ 

so that /?infl w o^. Thus, cuts in constant /?infl/t^res systemati- 
cally remove BH masses from the low-mass and low- velocity- 
dispersion portion of the M-a plane. For this reason, we ex- 
clude only those BH mass measurements that we believe to 
be unreliable. 

(7) Implications for space density of BHs. Our findings 
that the intrinsic scatter in the M-a for ellipticals and M- 
L relations for early-type galaxies areeo = 0.31±0.06 and 
eo = 0.38 ± 0.09, respectively, have an important influence 
on the determination of the space density of the most mas- 
sive BHs. We find that the density of BHs with Mbh > 
3 X 10^ Mg is 8 X 10"^ Mpc"^ (based on the M-a relation) 
and ~ 3 X 10~^ Mpc~^ (based on the M-L relation). 
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Institute of Technology, under contract with NASA. Finan- 
cial support was provided by NASAJHST grants GO-5999, 
GO-6587, GO-6633, GO-7468, and GO-9107 from the Space 
Telescope Science Institute, which is operated by AURA, Inc., 
under NASA contract NAS 5-26555. 



APPENDIX 
A. FITTING METHOD 

In this Appendix we describe our fitting method. Previous methods used for this problem include the Akritas & Bershady 
(1996) extension of the least-squares estimator (Ferrarese & Merritt 2000) and a symmetric method (Tremaine et al. 2002). 
Neither of these methods, however, can incorporate upper limits naturally. We use two general methods for fitting the M-a 
relation in this paper: (1) the symmetric %^ method of Tremaine et al. (2002), and (2) a generaUzed maximum-likelihood method 
developed for this work. We include the symmetric method for comparison to Tremaine et al. (2002), but the maximum 
likelihood method more naturally includes an intrinsic scatter component and more naturally incorporates upper limits. 

We consider several functional forms for the measurement error distribution and for the intrinsic scatter, both because it is far 
from clear that the conventional log-normal distribution accurately describes either distribution, and because distributions with 
fatter tails than normal tend to handle outUers more robustly. In §§ A.2-A.6 we describe our various assumptions about the error 
distribution and scatter. 

Ultimately, based on tests described in Appendix B below, we find that the choice of error distribution does not significantly 
change the fitted slope, intercept, or scatter; thus, we adopt the results from the fit with Gaussian distributions as the most 
straightforward. 

We take the measurement error in logarithmic mass to be 0.5[log(Mhigh) — log (Miovf)], where Mhigh andMiow are the pubUshed 
bounds to the la range in BH mass. Even though some of the galaxies in our sample have asymmetric errors (i.e., Mbh — Mow 
Mhigh — Mbh), we interpret the errors to be symmetric in logarithmic mass. This is a shortcoming of our methods, which should 
be remedied in future work. For the sake of consistency, we interpret the errors to be the 68% confidence intervals regardless of 
the form of the error distribution, as described below. In our sample, IC 1459 and NGC 1399 each have two mass measurements 
that are marginally inconsistent yet reliable. To account for this we include both measurements and weight each half as much. 

Similar to the measurement errors in BH mass, we take the magnitude of the intrinsic scatter to be the interval that contains 
68% of the area under the curve, regardless of the shape of the intrinsic scatter. We further assume (1) that the shape of the 
intrinsic scatter is synometric above and below the ridge Une, and (2) that the magnitude of the intrinsic scatter is independent of 
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Fig. 14. — BH space density derived from velocity dispersion (left) and luminosity (right) of the host galaxy. The curves show logarithmic cumulative density 
for different assumptions. The curves are the BH mass functions based on {left) the SDSS velocity dispersion function (Sheth et al. 2003) and (right) a fit of SDSS 
data to the Schechter (1976) luminosity function fit (Blanton et al. 2003) augmented by the brightest cluster galaxies of Postman & Lauer (1995). Predictions 
of BH mass based on velocity dispersion (left) and luminosity function (right) are from fits to all galaxies in this paper (orange) and to elliptical galaxies in this 
paper (red). We also plot corresponding predictions from Lauer et al. (2007b) (blue). Dashed lines indicate mass functions assuming no intrinsic scatter, and 
solid lines show mass functions assuming the best-fit value of the intrinsic scatter. For reference, both panels show alternate density calculations. The magenta 
dot-dashed line is based on the Richards et al. (2005) luminosity function from SDSS at a redshift of z = 2.5, assuming (1) a bolometric correction of 6.5 (Elvis 
et al. 1994; Sazonov et al. 2004), (2) that the AGN are accreting at their Eddington limit, and (3) that the AGN have a duty cycle of 0.03 (Steidel et al. 2002). The 
purple dot-dashed line is the quasar model due to Hopkins et al. (2006a). The differences between the solid and dashed curves of the same color show the large 
discrepancy between considering and omitting intrinsic scatter in estimating the density of the BHs. The red, solid curves in both panels show the best-fit results 
from this paper. 



velocity dispersion. 

A. 1 . Maximum-Likelihood Method for Parameter Estimation 
Let P be the probability that there is a BH of logarithmic mass /j = log (M/M . ) given galaxy properties s. For our fits, s is 
either the logarithm of the velocity dispersion s = log (o/200 km s^') or the logarithm of the luminosity s = log (Lv/10"L. y)- 
Thus, P{fj\s) is either the M-a or M-L relation. Given a set of observations chosen on the basis of s (or other properties not 
correlated with BH mass), the likelihood of a set of observed points {^u, ,*,} is the product of the likelihood of each pair of 
measurements ij — £{/Ji,Si): 

L^Ylli. (Al) 

In the absence of measurement errors 

li^P{lii\si). (A2) 

For a given observation, the probability of measuring a mass between /Jobs and /Job,s + '^A'obs in galaxy /, given that the actual BH 
mass in this galaxy is /7, is 2i(A'obs|A')<^A'obs, normalized so that the integral over all /Job,s is unity. Then equation (A2) becomes 

ii^eifii.st)^ f Qi{iii\^i)P{ji\si)dii. (A3) 

Jail Al 

Note that equation (A3) reduces to equation (A2) if = 5(/j — /j,). We also include observational upper limits to BH 

masses in our fits. In doing so, we must also allow for the possibility that some galaxies have no BH; otherwise, a single galaxy 
with no BH and a strong observational upper limit could strongly bias our parameter fits. To take this into account, we modify 
equation (A3) to be 

^, = P0G,(A'd + -^0) / Qi{^ii\^i)P{^i\s,)d^i, (A4) 

Jail Al 

where fa is the probability of the galaxy having no BH. In all cases is consistent with zero, with typical estimates of P% = 
0.003 ±0.03. 

In the following sections we consider specific forms for ^nd P(/j|s,). 
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TABLE 2 
List of Abbreviations Used 



— — 

Abbrev. 


Qi(p) 




— — 

Short Description 


GG 


A9 


A5 


Gaussian error distribution with Gaussian intrinsic scatter 


CG 


A12 


A5 


Constant probability errors with Gaussian intrinsic scatter 


DG 


A14 


A5 


Double-sided exponential errors with Gaussian intrinsic scatter 


DD 


A14 


A15 


Double- sided exponential errors 








with double-sided exponential intrinsic scatter 


LG 


A16 


A5 


Lorentzian errors with Gaussian intrinsic scatter 


LL 


A16 


A17 


Lorentzian errors with Lorentzian intrinsic scatter 


SU 






Full sample with upper limits 


S 






Sample without upper limits 


RS 






Restricted sample, no upper limits 



Note. — List of abbreviations. Column 1 gives the abbreviation. For the maximum- 
likelihood method abbreviations. Columns 2 and 3 give the equation numbers that describe the 
measurement error distribution Qiip), and the intrinsic scatter distribution Pi(p), respectively, 
while Colvmin 4 gives a shgrt description. 

A.2. Gaussian Error Distribution with Gaussian Scatter 



LetGe(xi— X2) = (l/e\/27t)exp[— (xi — X2)^/2e^] be a nonnalized Gaussian with dispersion e. Then suppose 

P{ti\si) = G,,{^-f{si)), (A5) 

where / is a ridge line through the joint distribution in /v and s. That is, / is the M-o or the M-L relation that we are considering. 
We refer to eo as the cosmic or intrinsic scatter in the relation, and we assume it to be independent of all galaxy properties. Then 



A = J Qdl^,\l^)G^,iiJ-f{s,))diJ. (A6) 

If Qi is normally distributed about jJi (i.e., the error in measured logarithmic BH mass is Gaussian) so that 

Qi = G,,{^^i-^i), (A7) 

where e,- is the measurement error, then 

d = I Ge. ini -n) Geo {M-f{si))dM. (A8) 
Equation (A8) is a convolution of two Gaussians — itself a Gaussian with variance equal to the sum of the two variances: 

^«- = G^/^(/i,- /(.,■))• (A9) 

Equation (A9) is the justification for adding cosmic scatter to the measurement error in Tremaine et al. (2002). We refer to this 
model as Gaussian error distribution with Gaussian intrinsic scatter (GG; see Table 2 for a hst of abbreviations for models and 
samples). This is the method we use for our final values of the parameters in the M-a and M-L relations. 

A.3. Upper Limits 

To include upper limits, we note that the probability that a galaxy with properties s has a logarithmic BH mass greater than 

is 

U{^,u\s) = {l-Pf>) rP{^\s)d^. (AlO) 

Suppose that an observation indicates the mass of the BH in a given galaxy is less than jU„ at the nagg level. This means that the 
mass is greater than /j„ with probability 5„ and less than jU„ with probability 1 — 5„, where 5i = 0.159, 82 = 0.0228, 83 = 0.00135, 
etc. To include the upper limits in the maximum-likelihood method, we must include the possibility that the observation has 
incorrectly concluded that the mass is less than which will happen a fraction 5„ of the time. Then the hkelihood of observing 
an upper limit of mass /i„ at the nagg level in a galaxy with property s is 

li = KU{nu\s) + (1 -5„)[1 - U{nu\s)]. (All) 
A.4. Constant Probability 

We may also decide that a mass measurement means the mass is restricted between ii\ and /Y2, but that there is no preferred 
mass in that range. We may assume this, for example, because we believe that the published error bars are a reasonable estimate 
of the uncertainty in BH mass, but the error distribution is unknown. In this case, Qi is a constant C, and 

li = C Ge,{fi-f{Si))dfi (A12) 

with 

C= — ^- — , (A13) 

where the normalization arises since the integral of Qi over all jU, should be unity. In order to maintain the consistency of having 
e,- indicate the 68% confidence Umit, we take /ii, 2 = W T e(/0.68. We refer to this method as CG. 
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Fig. 15. — Comparison of assumed error distributions with measured probability distribution of NGC 4026 as a function of Mbh (l^ff) and log (Mbh/Mt ) 
(right). The measured probability distribution comes from the chP' values from Giiltekin et al. (2009). The distributions are normalized so that the integral over 
log (MBH/M0) is unity. This figure illustrates the differences among the assumed error distributions. 

A.5. Double-Sided Exponential 

We also use robust methods to estimate the parameters of /(i,) by dropping the assumption of a Gaussian distribution of 
measurement errors in /j and, optionally, the assumption of Gaussian intrinsic scatter in /j. Robust methods are, in general, more 
tolerant of outliers in distributions. This may be especially appropriate for our sample since the measurements come from several 
different groups using several different measurement methods, and the systematic errors may be larger than the intrinsic scatter. 
There is no a priori reason to assume that the intrinsic scatter is correctly described by a Gaussian distribution. We use two robust 
methods. First, we use a double-sided exponential, a standard robust method (e.g.. Press et al. 1986): 

Qi{lii\^i) = (2fl,)"' exp(- \ij-ij,\/ai), (A14) 

where we determine a, by assuming that the quoted logg measurement errors contain the 68% confidence interval and thus 
correspond to — In (0.32)fl,. We use this form of measurement error with either equation (A5) (DG) or with a cosmic scatter that 
is also described by a double-sided exponential (DD): 

Pi{p) = {2ao)-^exp{-\pi-f{si)\/ao), (A15) 

where ao = —Eq/ In 0.32. 

A.6. Lorentzian 

For our second robust method we use a Lorentzian to describe the error distribution of the BH mass measurements. Lorentzian 
distributions have large tails and thus tend to be especially tolerant of outliers: 

where we determine F, by assuming the 68% confidence intervals correspond to tan(0.687t/2)F,. We use this form with either 
equation (A5) (LG) or with a cosmic scatter that is also described by a Lorentzian (LL): 

FoA l_ 



PM^\^]-, (A17) 



where Fq = eo/ tan (0.68:1/2). 

The error distributions described above are compared to the measured probability found for NGC 4026 by Giiltekin et al. (2009) 
in Figure 15. 

In practice, for all methods we minimize — IniL by the downhill simplex method (e.g.. Press et al. 1986). To calculate uncer- 
tainties in parameters, we change each parameter and refit (allowing the other parameters to vary to maximize InX) until IniL 
decreases by 0.5. We also run a Monte Carlo bootstrap described in Appendix B. 

A.7. Errors in the Independent Variable 

Errors in the independent variables (velocity dispersion and bulge luminosity) are incorporated by Monte Carlo sampling of 
the independent variables. For example, for determining the M-o fit, we select values of a according to the measured values and 
uncertainties. For each fit 10^ realizations are done, and the small parameter uncertainties from this are added in quadrature. 

A.8. Maximum-Likelihood Method: Model Comparison 
To compare two models we calculate the odds ratio 

_ J La{ai,a2,...,a,n)Pa{ai,a2,...,a,„)daida2...da,n (Al^) 
J Lb{bi,b2,-..,b„)Pi,{bi,b2,—,b„)dbidb2...db„ ' 
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TABLE 3 

Parameter Estimates for M-a Relation 



Method 


Sample 


a 


B 


EO 




T02 








41 




1 uzey 


O 

o 




a QQ _i_ '^/^Q 








c 


S 1 Q -1- ft C\f,A 


4 ft^i -1- ft ^7ft 












A 94 _|_ 4, 1 ft 


zlzl -1- OSQ 


OflO/l -1-0 018 


CG 


su 


8.13±0.085 


4.28 ±0.437 


0.45 ± 0.063 


0.0006 ±0.016 


DG 


su 


8.09 ±0.088 


4.37 ±0.603 


0.52 ±0.064 


0.0002 ±0.015 


LJLJ 




o . i o It U.U / J 


A -1- ft '?89 


ft Aft -1- ft ftf^Q 


ft ftftft 1 -1- n ft 1 '^ 


LG 


SU 


o . 1 J zn yj .KJ 1 " 




-1- ft ftfSS 


ft ftftft^ -1- ft? 1 


LL 






A ftft -1- ft AS^fy 


-1- ft 1 OS 


ft ftftft? -1- n ft?n 




c 
o 


J3 1 Q _|_ ft7Q 


a QC _i_ ATyX 


ft A'^ -1- ft ftSR 




CG 


s 


8 18 + 079 


3 96 d- 426 


v<^J —1— VJ-vJU 




DG 


s 


8. 15 ±0.093 


4.05 ±0.507 


0.51 ±0.067 




DD 


s 


8.23 ±0.073 


3.88 ±0.760 


0.39 ±0.082 




LG 


S 


8.21 ±0.073 


3.91 ±0.676 


0.37 ±0.068 




LL 


s 


8.27 ±0.072 


3.71 ±0.402 


0.32 ±0.094 




GG 


RS 


8.29 ±0.078 


3.74 ±0.404 


0.25 ±0.059 




CG 


RS 


8.30 ±0.069 


3.76 ±0.369 


0.25 ±0.050 




DG 


RS 


8.29 ±0.073 


3.71 ±0.405 


0.24 ±0.059 




DD 


RS 


8.33 ±0.067 


3.73 ±0.342 


0.20 ±0.060 




LG 


RS 


8.30 ±0.083 


3.72±0.417 


0.23 ±0.070 




LL 


RS 


8.38 ±0.087 


3.74 ±0.708 


0.15±0.108 





Note. — Results from fits. The first three lines are the methods of Tremaine et al. 
(2002). The first line (T02) is the average of two symmetric variants in which intrinsic 
scatter is increased until per degree of freedom is unity, with each mass measurement 
either given equal weight (T02eq) or with individual errors (T02ind). Method and sam- 
ple abbreviations are described in Table 2. 



where La is the hkehhood of the data given model a with m parameters a,-, which have a prior probabihty distribution Pa, 

and similarly for model b. The models need not have the same number of parameters. We assume that the prior probability 
distributions are uniform within the ranges [4, 12] for intercepts \ju{pe = 200 km s"') ox ^l{Lv = IO^^Lq^v)], [-10, 10] for slopes, 
and [0, 2] for cosmic scatter, though any reasonable set of ranges produces the same qualitative results. The integrals are calculated 
by the Vegas Monte Carlo method (Lepage 1978). 

B. RESULTS FROM DIFFERENT ERROR DISTRIBUTIONS, SCATTER FORMS, AND FIT SAMPLES 

The results from M-a fits are summarized in Table 3 and plotted in Figure 16. 

First, we analyze our dataset excluding upper limits (S) with the method of Tremaine et al. (2002), which consists of two distinct 
variants. (1) Measurement errors and intrinsic scatter are combined into a single number, the rms deviation of log(MBH/M0) 
from the ridge line (denoted eo, even though this is normally reserved for the intrinsic scatter), which is assumed to be the same 
for all galaxies and is determined by requiring that the reduced of the fit be unity. (2) The stated measurement errors and 
the intrinsic scatter eo are added in quadrature, as in Equation (A9). When using the first method, we find a = 8. 19 ±0.06 and 
P = 3.99 ± 0.37 for 8o = 0.43. When using the second method, we find a = 8. 19 ± 0.06 and p = 4.06 ± 0.37 for eo = 0.40. The 
slopes we obtain here are consistent with Tremaine et al. (2002), but the scatters are larger (see § 5.2). 

Then we use the maximum-likehhood method with the various assumed forms for measurement errors and the intrinsic scatter. 
Using the different error distributions for the sample that includes upper limits (SU), we find a range of intercepts from 8.09 to 
8.23, which is consistent compared to the largest parameter uncertainty estimate of 0.09. Similarly, we find slopes in the range 
4.00 to 4.37, a smaller range than the ~ 0.4 uncertainty. The estimates of the intrinsic scatter range from 0.35 to 0.52 with 
uncertainties of 0.06-0.1 1. 

For the sample that does not include upper limits (S), we find a similar self-consistency among the different error and scatter 
distributions. The intercept is larger when excluding upper limits. The slope is shallower in the sample without limits, though not 
at a significant level, and the scatter is slightly (but not significantly) smaller. The consistency of the results among the different 
error distributions for a given sample indicates that the choice of distribution is not driving the values we obtain for parameters. 
The results for the restricted sample (RS) are similarly self-consistent, though compared to the full sample SU they have a larger 
intercept, shallower slope, and smaller intrinsic scatter, for reasons discussed in § 4. Thus, assuming a different error distribution 
does not have a strong effect on the results, but selecting a different sample does. 

As a check on the uncertainty estimation, we calculate a bootstrap Monte Carlo of our sample. For each combination of 
measurement error and intrinsic scatter distributions, we randomly extract data points from the original sample with replacement 
for 10^ reahzations and fit each sample. The bootstrap method provides a powerful way of examining whether any outlying 
measurements are driving the fit as well as a consistency check on uncertainty estimations. We show the distribution of best- fit 
parameters for two combinations of distributions: Gaussian measurement error with Gaussian intrinsic scatter (Fig. 17) and the 
Lorentzian measurement error with Lorentzian intrinsic scatter (Fig. 18), each with and without upper limits. Table 4 shows that 
in all cases the median fit parameters with 68% intervals of the distributions are very close to the values we found above. The 
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Intercept Slope 



Fig. 16. — Results of fits to the M-a relation showing best-fit intercepts (a), best-fit slopes (p), and best-fit intrinsic scatters (Eq) against each other. The 
adopted fitting methods are plotted with the symbol from the legend, which uses abbreviations from Table 2. For clarity, we show representative error ellipses for 
method GG for each of the samples. The error ellipses are the 68% confidence limits in the joint distribution of the two parameters plotted. The sample used is 
indicated by the color of the parameter uncertainties: blue (sample without limits [S]), red (full sample including upper limits [SU]), and green (restricted sample 
[RS]). The parameter values for the various adopted methods do not significantly vary, as can be seen from the fact that most points from the same sample (same 
color) fall inside the uncertainty box implied by the error bars. The values from different samples, however, do vary. We discuss the reasons for the differences 
in sample RS (which has a cutoff in Rini/dres) in § 4 and for the differences when including upper limits in § 5.4. 



TABLE 4 
Bootstrap Monte Carlo 



Method Sample a 



GG SU 

LL SU 

GG S 

LL S 



P EO 



8.14±.07 4.31 ±.52 0.44±.07 

8.23 ±.08 3.96 ±.48 0.33 ±.09 

8. 17 ±.07 3.95 ±.36 0.42 ±.07 

8.27 ±.08 3. 78 ±.52 0.31 ±.08 



Note. — Results from bootstrap fits expressed as median 
value. The uncertainties quoted encompass 68% of the val- 
ues. Methods are as given in Table 3. The values are in very 
close agreement with the corresponding values given in Ta- 
ble 3. This indicates that the uncertainties deiived from our 
fitting method are accurate and are not strongly influenced by 
outliers. 



distributions appear well approximated by a Gaussian. 

Given the consistency among the different assumptions and the consistency between the different uncertainty estimation meth- 
ods, we conclude that the results are insensitive to the choice of distribution of measurement error and intrinsic scatter (so long 
as the "width" Eq is consistently defined as the interval containing 68% of the distribution). The best-fit parameters of the M-a 
and M-L relations, however, do depend on the sample choice. Thus, we adopt the simplest combination of measurement error 
and intrinsic scatter distributions, GG, for the estimates of the fit parameters and their uncertainties that we quote in the abstract. 



Intrinsic Scatter of the M-a and M-L Relations 
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Fig. 17. — Histograms of parameter estimates from the bootstrap samples assuming Gaussian measurement eiTors and Gaussian intrinsic scatter. Left panels 
are for the full sample, SU; right panels do not include upper limits (sample S). The distributions are unimodal, and the medians and 68% intervals are in good 
agreement with our best-fit values and uncertainties. 




Fig. 18. — Histograms of parameter estimates from the bootstrap samples assuming Lorentzian measurement en'ors and Lorentzian intrinsic scatter. Left panels 
are for the full sample, SU; right panels do not include upper limits (sample S). The distributions are unimodal, and the medians and 68% intervals are in good 
agreement with our best-fit values and uncertainties. 



22 



Giiltekin et al. 



REFERENCES 



Adams, R C, Graff, D. S., Mbonye, M., & Richstone, D. O. 2003, ApJ, 591, 
125 

Adams, R C, Graft', D. S., & Richstone, D. O. 2001, ApJ, 551, L31 
Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 
Atkinson, J. W., et al. 2005, MNRAS, 359. 504 

Barth, A. J., Sarzi, M., Rix, H.-W., Ho. L. C, Filippenko, A. V.. & Sargent, 

W. L. W. 2001, ApJ, 555, 685 
Beifiori, A., Sarzi, M., Corsini, E. M., Dalla Bonta, E., Pizzella, A., Coccato, 

L., & Bertola, R 2008, preprint (0809.5103) 
Bender, R., et al. 2005, ApJ, 631, 280 
Bentz. M. C, et al. 2006, ApJ, 651, 775 
Bernardi, M., et al. 2006, AJ, 131, 2018 
Blanton, M. R., et al. 2003, ApJ, 592, 819 
Bower, G. A., et al. 1998, ApJ, 492, LI 1 1 

. 2001, ApJ, 550, 75 

Burkert, A., & Silk, J. 2001, ApJ, 554, L151 
Burstein, D. 1979, ApJ, 234, 435 

Capaccioli, M., Held, E. V., & Nieto, J.-L. 1987, AJ, 94, 1519 

Capetti, A., Marconi, A., Macchetto, D., & Axon, D. 2005, A&A, 431, 465 

Cappellari, M., Neumayer, N., Reunanen, J., van der Werf, R R, de Zeeuw, 

R T., & Rix, H. . 2008, preprint (0812.1000) 
Cappellari, M., Verolme, E. K., van der Marel, R. P., Kleijn, G. A. V., 

lUingworth, G. D.. Franx. M., CaroUo, C. M., & de Zeeuw, R T. 2002, 

ApJ, 578, 787 

Coccato, L., Sarzi, M., Pizzella, A., Corsini, E. M., Dalla Bonta, E., & 

Bertola, R 2006, MNRAS, 366, 1050 
Cretton, N., & van den Bosch, R C. 1999, ApJ, 514, 704 
Dalla Bonta, E., Fen'arese, L., Corsini, E. M., Miralda-Escude, J., Coccato, 

L., Sarzi, M., Pizzella, A., & Beifiori, A. 2008, preprint (0809.0766) 
de Francesco, G., Capetti, A., & Marconi, A. 2006. A&A, 460, 439 
. 2008, preprint (0801.0064) 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Jr., Buta, R. J., 
Paturel, G., & Fouque, P. 1991, Third Reference Catalogue of Bright 
Galaxies (Springer- Verlag Berlin Heidelberg New York) 

Devereux, N., Ford, H., Tsvetanov, Z., & Jacoby, G. 2003, AJ, 125, 1226 

Dressier, A. 1989, in lAU Symposium 134, Active Galactic Nuclei, ed. D. E. 
Osterbrock & J. S. Miller, 217 

Elvis, M., et al. 1994, ApJS, 95, 1 

Emsellem, E., Dejonghe, H., & Bacon, R. 1999, MNRAS, 303, 495 
Erwin, R 2004, A&A, 415, 941 

Fen'arese, L. 2002, in Current high-energy emission around black holes, ed. 

C.-H. Lee & H.-Y. Chang, 3 
Ferrarese, L., & Ford, H. C. 1999, ApJ, 515, 583 
Ferrarese, L., Ford. H. C, & Jaffe, W. 1996, ApJ, 470, 444 
Fen'arese, L., & Ford, H. 2005, Space Science Reviews, 116, 523 
Fenarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Gebhardt, K., & Thomas, J. 2009, ApJ, submitted, 
Gebhardt. K., et al. 2000a, ApJ, 539, L13 

. 2000b, ApJ, 543, L5 

. 2000c, AJ, 119, 1157 

. 2001, AJ, 122, 2469 

. 2003, ApJ, 583, 92 

. 2007, ApJ, 671, 1321 

Ghez, A. M., Salim, S., Homstein, S. D., Tanner, A., Lu, J. R., Monis, M., 

Becklin, E. E., & Duchene, G. 2005, ApJ, 620, 744 
Ghez, A. M., et al. 2008, ApJ, 689, 1044 

Gillessen, S., Eisenhauer, R, Trippe, S., Alexander, T, Genzel, R., Martins, 

R, & Ott, T. 2008, preprint (0810.4674) 
Gnedin, O. Y, Kravtsov, A. V., Klypin, A. A., & Nagai, D. 2004, ApJ, 616, 

16 

Graham, A. W. 2008, ApJ, 680, 143 

Graham, A. W., Erwin, P, Caon, N., & Trujillo, I. 2001, ApJ, 563, LU 
Greenhill, L. J., & Gwinn, C. R. 1997, Ap&SS, 248, 261 
Greenhill, L. J., Moran, J. M., & Herrnstein, J. R. 1997, ApJ, 481, L23 
Greenhill, L. J., et al. 2003, ApJ, 590, 162 
Giiltekin, K., et al. 2009, ApJ, accepted, 

Haring-Neumayer, N., Cappellari, M., Rix, H.-W., Hartung, M., Prieto, 
M. A., Meisenheimer, K., & Lenzen, R. 2006, ApJ, 643, 226 

Herrnstein, J. R., Moran. J. M., Greenhill, L. J., & Trotter, A. S. 2005, ApJ, 
629,719 

Ho, L. C, Sarzi, M., Rix, H.-W., Shields, J. C, Rudnick, G., Pilippenko, 

A. v., & Barth, A. J. 2002, PASP, 1 14, 137 
Hopkins, P. F, Hernquist, L., Cox, T. J., Robertson, B., & Springel, V. 2006a, 

ApJS, 163, 50 

Hopkins, P. R, Robertson, B., Krause, E., Hernquist, L., & Cox, T. J. 2006b, 
ApJ, 652, 107 



Houghton, R. C. W., MagoiTian, J., Sarzi, M., Thatte, N., Davies, R. L., & 

Krajnovic, D. 2006, MNRAS, 367, 2 
Hu, J. 2008, MNRAS, 386, 2242 

Knapen, J. H., Mazzuca, L. M., Boker, T, Shlosman, I., Colina, L., Combes, 

R, & Axon, D. J. 2006, A&A, 448, 489 
Kormendy, J. 1988, ApJ. 335, 40 

. 1993, in The Nearest Active Galaxies, ed. J. Beckman, L. Colina, & 

H. Netzer (Madrid: Consejo Superior de Investicaciones Cientificas), 197 

. 2004, in Coevolution of Black Holes and Galaxies, ed. L. C. Ho 

(Cambridge: Carnbridge Univ. Press), 1 

Kormendy, J., Pisher, D. B., Cornell, M. E., & Bender, R. 2008, preprint 
(0810.1681) 

Kormendy, J., & Gebhardt, K. 2001, in American Institute of Physics 
Conference Series 586, 20th Texas Symposium on relativistic astrophysics, 
ed. J. C. Wheeler & H. Martel (Melville, NY: AIP), 363 

Kormendy, J., & Illingworth, G. 1982, ApJ, 256, 460 

. 1983, ApJ, 265, 632 

Kormendy, J., & Kennicutt. R. C, Jr. 2004, ARA&A, 42, 603 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 
Kormendy, J., et al. 1997, ApJ, 482, L139 
Krajnovic, D., Sharp, R., & Thatte, N. 2007, MNRAS, 374, 385 
Lauer, T. R., Trernaine, S., Richstone, D., & Raber, S. M. 2007a, ApJ, 670, 
249 

Lauer, T. R., et al. 1996, ApJ, 471, L79 

. 2005, AJ, 129, 2138 

. 2007b, ApJ, 662, 808 

. 2007c, ApJ, 664, 226 

Lepage, G. P. 1978, Journal of Computational Physics, 27, 192 
Lodato, G., & Berlin, G. 2003, A&A, 398, 517 

Macchetto, P., Marconi, A., Axon, D. J., Capetti, A., Sparks, W., & Crane, P. 

1997. ApJ, 489, 579 
Maciejewski, W., & Binney, J. 2001, MNRAS, 323, 831 
Magonian, J. 2006, MNRAS, 373, 425 
Magonian, J., et al. 1998, AJ, 115, 2285 

Marconi, A., Capetti, A., Axon, D. J., Koekemoer, A., Macchetto, D., & 

Schreier, E. J. 2001, ApJ, 549, 915 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 

Marconi, A., Pastorini, G., Pacini, P., Axon, D. J., Capetti, A., Macchetto, D., 
Koekernoer, A. M., & Schreier, E. J. 2006, A&A, 448, 921 

Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R., & Salvati, M. 
2004, MNRAS, 351, 169 

Marconi, A., et al. 2003, ApJ, 586, 868 

Martel, A. R., Turner, N. J., Sparks, W. B., & Baum, S. A. 2000, ApJS, 130, 
267 

Merritt, D., & Ferrarese, L. 2001a, in Astronomical Society of the Pacific 
Conference Series 249, The Central Kiloparsec of Starbursts and AGN: 
The La Raima Connection, ed. J. H. Knapen, J. E. Beckman, L Shlosman, 
& T. J. Mahoney, 335 

. 2001b, ApJ, 547, 140 

Menitt, D., Ferrarese, L., & Joseph, C. L. 2001, Science, 293, 1 1 16 
Miyoshi, M., Moran, J., Hennstein, J., Greenhill, L., Nakai, N., Diamond, P., 

& Inoue, M. 1995, Nature, 373, 127 
Neumayer, N., Cappellari, M., Reunanen, J., Rix, H.-W., van der Werf, P. P., 

de Zeeuw, R T, & Davies, R. L 2007, ApJ, 671, 1329 
Novak, G. S., Paber, S. M., & Dekel, A. 2006, ApJ, 637, 96 
Nowak, N., Saglia, R. P, Thomas, J.. Bender, R., Pannella, M., Gebhardt, K., 

& Davies, R. I. 2007, MNRAS, 379, 909 
Onken, C. A., Fenarese, L., Meiritt, D., Peterson, B. M., Pogge, R. W., 

Vestergaard, M., & Wandel, A. 2004, ApJ, 615, 645 
Onken, C. A., et al. 2007, ApJ, 670, 105 
Pastorini, G., et al. 2007, A&A, 469, 405 

Peng, C. Y, Impey, C. D., Rix, H.-W., Kochanek, C. S., Keeton, C. R., Falco, 

E. E., Lehar, J., & McLeod, B. A. 2006, ApJ, 649, 616 
Peterson, B. M., et al. 2004, ApJ, 613, 682 
Postman, M., & Lauer, T. R. 1995, ApJ, 440, 28 

Press, W. H., Flannery, B. P., & Teukolsky, S. A. 1986, Numerical Recipes. 

The Art of Scientific Computing (Cambridge: Cambridge University Press) 
Richards, G. T, et al. 2005, MNRAS, 360, 839 

Richstone, D. 2004, in Coevolution of Black Holes and Galaxies, ed. L. C. 

Ho (Cambridge: Univ. Chicago Press), 280 
Richstone, D., et al. 1998, Nature, 395, A14 

Sarzi, M., Rix, H.-W., Shields, J. C, Rudnick, G., Ho, L. C, Mcintosh, D. H., 

Filippenko, A. V., & Sargent. W. L. W. 2001, ApJ, 550, 65 
Sarzi, M., et al. 2002, ApJ, 567, 237 

Sazonov, S. Y, Ostriker, J. R, & Sunyaev, R. A. 2004, MNRAS, 347, 144 
Schechter, R 1976, ApJ, 203, 297 



Intrinsic Scatter of the M-a and M-L Relations 



23 



Shen, J., Vanden Berk, D. E., Schneider, D. P., & Hall, P. B. 2008, AJ, 135, 
928 

Shen, Y., et al. 2007, AJ, 133, 2222 
Sheth, R. K., et al. 2003, ApJ, 594, 225 

Silge, J. D., Gebhardt, K., Bergmann, M., & Richstone, D. 2005, AJ, 130, 
406 

Silk, J., & Rees, M. J. 1998, A&A, 331, LI 
Siopis, C, et al. 2008, preprint (0808.4001) 

Steidel, C. C, Hunt, M. R, Shapley, A. E., Adelberger, K. L., Pettini, M., 

Dickinson, M., & Giavalisco, M. 2002, ApJ, 576, 653 
Stephens, M. 1974, Journal of the American Statistical Association, 69, 730 
Tadhunter, C, Marconi, A., Axon, D., Wills, K., Robinson, T. G., & Jackson, 

N. 2003, MNRAS, 342, 861 
Tonry, J. L., Dressier, A., Blakeslee, J. P., Ajhar, E. A., Fletcher, A. B., 

Luppino, G. A., Metzger, M. R., & Moore, C. B. 2001, ApJ, 546, 681 
Tremaine, S. 1995, AJ, 1 10, 628 
Tremaine, S., et al. 2002, ApJ, 574, 740 



Treu, T., Malkan, M. A., & Blandford, R. D. 2004, ApJ, 615, L97 
Treu, T., Woo, J.-H., Malkan, M. A., & Blandford, R. D. 2007, ApJ, 667, 117 
Valluri, M., FeiTarese, L., Merritt, D., & Joseph, C. L. 2005, ApJ, 628, 137 
Valluri, M., Merritt, D., & Emsellem, E. 2004, ApJ, 602, 66 
van den Bosch, F. C, Jaffe, W., & van der Marel, R. P 1998, MNRAS, 293, 
343 

van der Marel, R. P, & van den Bosch, F. C. 1998, AJ, 1 16, 2220 
Verolme, E. K., et al. 2002, MNRAS, 335, 517 

Vestergaard, M., Fan, X., Tremonti, C. A., Osmer, P. S., & Richards, G. T. 

2008, ApJ, 674, LI 
Volonteri, M. 2007, ApJ, 663, L5 

Wold, M., Lacy, M., Kaufl, H. U., & Siebenmorgen, R. 2006, A&A, 460, 449 
Wyithe, J. S. B. 2006a, MNRAS, 365, 1082 

. 2006b, MNRAS, 371, 1536 

Yu, Q., & Tremaine, S. 2002, MNRAS, 335, 965 



24 Giiltekin et al. 



TABLE 5 

Sample of Dynamically Detected Black Hole Masses 
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3.4 X 


10* 


2.8 X 


10* 


4.9 X 


10* 


stars, 17 


213±10 


-21.88 


-21.80±0.20 


6.69 


RS 


N3607 


El 


19.9 


1.2 X 


10* 


7.9 X 


10'' 


1.6 X 


10* 


stars, 17 


229±11 


-21.62 


-21.62±0.10 


2.12 


RS 


N3608'' 


El 


23.0 


2.1 X 


10* 


1.4 X 


10* 


3.2 X 


10* 


stars, 6 


182±9 


-21.05 


-21.05±0.10 


2.17 


RS 


N3998 


SO 


14.9 


2.4 X 


10* 


6.2 X 


10' 


4.5 X 


10* 


gas, 18 


305±15 


-20.32 




1.52 


S 


N4026 


SO 


15.6 


2.1 X 


10* 


1.7 X 


10* 


2.8 X 


10* 


stars, 17 


180±9 


-20.28 


-19.83 ±0.20 


7.72 


RS 


N4258 


SABbc 


7.2 


3.78 X 


10'' 


3.77 X 


10'' 


3.79 X 


10' 


maser, 19 


115±10 


-21.31 




64.5 


RS 


N4261 


E2 


33.4 


5.5 X 


10* 


4.3 X 


10* 


6.6 X 


10* 


gas, 20 


315±15 


-22.72 


-22.72 ±0.06 


0.76 


RS 


N4291'' 


E2 


25.0 


3.2 X 


10* 


8.3 X 


10'' 


4.1 X 


10* 


stars, 6 


242±12 


-20.67 


-20.67±0.13 


1.41 


RS 


N43428 


SO 


18.0 


3.6 X 


10* 


2.4 X 


10* 


5.6 X 


10* 


stars, 21 


225±11 


-18.84 




0.35 


RS 


N4374I M84 


El 


17.0 


1.5 X 


lO' 


9.0 X 


10* 


2.6 X 


10' 


gas, 22 


296±14 


-22.45 


-22.45 ±0.05 


9.89 


S 


N4459'' 


E2 


17.0 


7.4 X 


10' 


6.0 X 


10' 


8.8 X 


10' 


gas, 12 


167±8 


-21.06 


-21. 06 ±0.04 


1.34 


S 


N4473'' 


E4 


17.0 


1.3 X 


10* 


3.6 X 


10' 


1.8 X 


10* 


stars, 6 


190±9 


-21.14 


-21.14±0.04 


2.11 


RS 


N4486 M87 


El 


17.0 


3.6 X 


lO' 


2.6 X 


10" 


4.6 X 


10' 


gas, 23 


375±18 


-22.92 


-22.92 ±0.04 


6.43 


RS 


N4486A 


E2 


17.0 


1.3 X 


10' 


9.0 X 


10* 


1.8 X 


10' 


stars, 24 


111±5 


-18.70 


-18.70±0.05 


3.74 


S 


N4564'' 


SO 


17.0 


6.9 X 


10' 


5.9 X 


10' 


7.3 X 


10' 


stars, 6 


162±8 


-20.10 


-19.60±0.32 


1.46 


RS 


N4594 


Sa 


10.3 


5.7 X 


10* 


1.7 X 


10* 


1.1 X 


10' 


stars, 25 


240±12 


-22.52 


-22.44 ±0.15 


8.48 


S 


N4596'' 


SBO 


18.0 


8.4 X 


10' 


5.9 X 


10' 


1.2 X 


10* 


gas, 12 


136±6 


-20.70 




1.86 


RS 


N4649'' M60 


E2 


16.5 


2.1 X 


10' 


1.5 X 


10« 


2.6 X 


10' 


stars, 6 


385±19 


-22.65 


-22.65 ±0.05 


10.2 


RS 


N4697'^ 


E6 


12.4 


2.0 X 


10* 


1.8 X 


10* 


2.2 X 


10* 


stars, 6 


177±8 


-21.29 


-21.29±0.11 


4.63 


RS 


N5077 


E3 


44.9 


8.0 X 


10* 


4.7 X 


10* 


1.3 X 


10' 


gas, 26 


222±11 


-22.04 


-22.04 ±0.13 


2.45 


RS 


N5128''" 


SO/E 


4.4 


3.0 X 


10* 


2.8 X 


10* 


3.4 X 


10* 


stars, 27 


150±7 


-21.82 


-21.82±0.08 


42.4 


S 


N5128j'" 


SO/E 


4.4 


7.0 X 


10^ 


3.2 X 


10' 


8.3 X 


10' 


stars, 32 


150±7 


-21.82 


-21.82±0.08 


42.4 


S 


N5576 


E3 


27.1 


1.8 X 


10* 


1.4 X 


10* 


2.1 X 


10* 


stars, 17 


183±9 


-21.26 


-21.26±0.13 


3.55 


S 


N5845* 


E3 


28.7 


2.9 X 


10* 


1.2 X 


10* 


3.4 X 


10* 


stars, 6 


234±11 


-19.77 


-19.77±0.13 


1.65 


S 


N6251 


El 


106.0 


6.0 X 


10* 


4.0 X 


10* 


8.0 X 


10* 


gas, 28 


290±14 






0.52 


S 


N7052 


E3 


70.9 


4.0 X 


10* 


2.4 X 


10* 


6.8 X 


10* 


gas, 29 


266±13 






0.65 


S 


N7457'' 


SO 


14.0 


4.1 X 


lO** 


2.4 X 


10* 


5.3 X 


10* 


stars, 6 


67±3 


-19.80 


-18.72±0.11 


0.53 


S 


N7582''8 


SBab 


22.3 


5.5 X 


10' 


4.4 X 


10' 


7.1 X 


10' 


gas, 30 


156±19'' 


-21.51 




0.22 


S 


A1836-BCG'' 


E 


157.5 


3.9 X 


10" 


3.3 X 


109 


4.3 X 


10' 


gas, 31 


288±14<' 


-23.31 


-23.31±0.15 


2.63 


S 


A3565-BCG'' 


E 


54.4 


5.2 X 


10* 


4.4 X 


10* 


6.0 X 


10* 


gas, 31 


322±16'' 


-23.27 


-23.27 ±0.15 


0.81 


S 
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Dist. 








Method, Oe 








Galaxy 


Type" 


Mpc 


Mo 


Me 


Ms 


Ref. km s"^ 




luO b 
'"V,bulge 


Rmn/dres Samp. 



References. — (1) Greenhill et al. (2003), (2) Cappellari et al. (2002), (3) Ghez et al. (2008) and Gillessen et al. (2008), (4) Verolme et al. (2002), (5) Bender et al. (2005), 
(6) Gebhardt et al. (2003), (7) Bower et al. (2001), (8) Lodato & Bertin (2003), (9) Atkinson et al. (2005), (10) Gebhardt et al. (2007), (1 1) Houghton et al. (2006), (12) Sarzi et al. 
(2001), (13) Devereux et al. (2003), (14) Emsellem et al. (1999), (15) Barth et al. (2001), (16) Gebhardt et al. (2000c), (17) Gultekin et al. (2009), (18) de Francesco et al. (2006), 
(19) Hermstein et al. (2005), (20) Ferrat^se et al. (1996), (21) Cretton & van den Bosch (1999), (22) Bower et al. (1998), (23) Macchetto et al. (1997), (24) Nowak et al. (2007), (25) 
Kormendy (1988), (26) de Francesco et al. (2008), (27) Silge et al. (2005), (28) Ferrarese & Ford (1999), (29) van der Marel & van den Bosch (1998), (30) Wold et al. (2006), (31) 
Dalla Bonta et al. (2008), (32) Cappellari et al. (2008). 

Note. — Mi„„ and My are the lower and upper limits of the allowed ranges in the BH mass measurement at the lo^g level. "Method" column indicates method of BH mass 
determination: stellar dynamics ("stars"), gas dynamics ("gas"), or maser dynamics ("maser"). Effective velocity dispersion is given as defined by equation (1). Dereddened 
magnitudes are given for the entire galaxy (AfJ j-) and the bulge (^^^l,|,lge)• The final column gives the most restricted sample of which each galaxy is a member. All members of the 
restricted sample (RS) are members of the sample without Umits (S), which are all members of the full sample with upper limits (SU). Galaxies that contain pseudobulges are marked 
by the superscript "g". 

a b 
Galaxy types are taken from de Vaucouleurs et al. (1991) with the following exceptions: NGC 3607, NGC 4459, and NGC 4564, which all come from Kormendy et al. (2008). Errors 

in bulge-disk decomposition are estimated from the range of decomposition values given in the literature and are propagated to an error in bulge luminosity. Bulge-disk decomposition 

is taken from Lauer et al. (2007b) with the following exceptions: NGC 1023 (Kormendy & Ilhngworth 1983), NGC 3115 (CapaccioU et al. 1987), and NGC 3384 (Burstein 1979).*^ The 

Circinus galaxy is in the plane of the Milky Way and has a BH mass measurement from masers (Greenhill et al. 2003). Ferrarese & Ford (2005) Ust the BH mass of Circinus as 

possibly in error because the incUnation of the maser disk is not constrained, yet the velocities obtained are clearly Keplerian, and the disk is unlikely to be far from edge-on. If 

the location of the maser detections indicates the extent of the disk, as has been assumed for NGC 4258 (Miyoshi et al. 1995), then the incUnation of the disk is probably ~ 75°, 

which woitld increase the mass by about 4% from the value Greenhill et al. (2003) find assixming edge-on inclination. We conclude that the presence of a BH is almost certain and 

that the mass and uncertainties obtained by Greenhill et al. (2003) are reliable estimates. If we omit Circinus from our sample, our best-fit parameter estimates are a — 8. 17 ±0.07, 

d e 
p — 4. 12 ±0.37, £o — 0.38 ±0.06. Central dispersion is given rather than effective 0^; see eq. (1). We use only the stellar dynamical measurement of the black hole in IC 1459. 

f 

The gas dynamical measurement may be in error because the gas kinematics are disturbed (M. Cappellari, private communication). The uncertainty in the mass of the Galaxy's black 
hole includes uncertainty in the distance to the Galactic center.^ The galaxy has a pseudobulge. If the classification is not obvious from the images available in NED or from the high 
bulge-to-disk ratio, we cite the reference for pseudobulge classification here: NGC 224 (Kormendy & lUingworth 1982), NGC 1300 (Knapen et al. 2006), NGC 2748 (Kormendy 
& Kennicutt 2004), NGC 2787 (Erwin 2004), NGC 3384 (Erwin 2004), NGC 4342 (van den Bosch et al. 1998).** Mass is 9% larger than originally pubUshed due to a numerical 
error.^ The BH in NGC 1 068 was originally measured to have Men = 1 5 x 1 O^M . by Greenhill & Gwinn (1997) under the assumption of Keplerian rotation in the disk despite 
the fact that they find that velocities fall off more slowly than Keplerian. Using the same data, Lodato & Bertin (2003) find a formally better fit assuming a massive gas disk, and, thus, 
we use their smaller value for Mbh. If we use the larger value, we obtain a = 8.13 ± 0.08, p = 4.21 ± 0.40, and eo = 0.44± 0.06.-' Both NGC 1399 and NGC 5128 have two different 
mass measurements that are at least mareinallv inconsistent with each other but appear to be individixally reliable. Rather than averaging these values we use both and weight each 
measurement half as much, as explained in Appendix A. The galaxies NGC 2787, NGC 4459, and NGC 4596 have masses reported by Sarzi et al. (2001) for both unconstrained 
disk inclination and for the best-fit disk inclination. We use the latter as the regions of the gas disk probed by the observations are likely to be aligned in the plane of the axisymmetric 
bulge (e.g. Mattel et al. 2000). These galaxies are all low velocity dispersion galaxies and thus unlikely to be triaxial.^ Maciejewski & Binney (2001) interpreted the double-peaked 
velocity profile in NGC 4374 as coming from a single component rather than from two separate components as Bower et al. (1998) did. The velocity profile is clearly doubly peaked 
as can be seen in figure 3 of Bower et al. (1998), and the second peak is likely to come from a more slowly rotating separate component, thus we include it in our sample. If the 
interpretation of Maciejewski & Binney (2001) is correct, the mass is a factor of 4 smaller than foimd by Bower et al. (1998). Omitting this galaxy, our best fit parameters change to 
a = 8. 12 ±0.08, P = 4. 19 ±0.41, andEo = 0.45 ±0.06.™ NGC 5128 has multiple Mbh measurements with ionized gas (Marconi et al. 200 1 ; Haring-Neumayer et al. 2006; Marconi 
et al. 2006), two-dimensional neutral gas velocity measurements (Krajnovic et al. 2007; Neumayer et al. 2007), and stellar dynamical measurements (Silge et al. 2005; Cappellari 
et al. 2008). Neumayer et al. (2007) argue that the ionized-gas measurements may be contaminated by the galaxy's jet. Three of the remaining measurements (Krajnovii et al. 2007; 
Neumayer et al. 2007; Cappellari et al. 2008) are in good agreement with each other, which we present as one value in addition to the other measurement (Silge et al. 2005). 
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TABLE 6 

Upper limits to black hole masses 



Dist. M„ Method, a, 



Galaxy 


Type 


Mpc 


Mo 


Confidence 


Ref. 




km s ^ 




1*0 a 

™F,bulge 




«mfi/"res 


Sample 


N3310 


SB(r)bc 


17.4 


4.2 X 10' 


2068 


gas. 


1 


83±4 


-20.56 






2.39 


SU 


N3351'' 


SBb 


8.7 


8.6 X 10* 


l068 


gas. 


2 


93±4 


-20.15 






0.90 


SU 


N3368 


SBab 


11.0 


3.7 X 10^ 


lOgg 


gas. 


2 


1I4±5 


-21.19 






1.80 


SU 


N3982 


SBb: 


18.2 


8.0 X 10' 


lags 


gas. 


2 


78±3 








8.75 


SU 


N3992 


SBbc 


18.2 


5.7 X 10' 


lags 


gas. 


2 


119±5 


-21.73 






1.95 


SU 


N4041 


S(rs)bc 


20.9 


6.4 X 10* 


3a6s 


gas. 


4 


88±4 


-20.40 






0.35 


SU 


N4143 


SBO 


16.8 


1.4 X 10* 


lags 


gas. 


2 


271±13 








1.59 


SU 


N4203 


SBO 


16.0 


3.8 X 10' 


loes 


gas. 


2 


110±5 


-20.34 






0.80 


SU 


N4321'' 


SBbc 


18.0 


2.7 X 10' 


l068 


gas. 


2 


74±3 


-21.95 






1.79 


SU 


N4435 


SBO 


17.0 


8.0 X 10* 


3068 


stars, 


5 


150±7 


-20.45 






0.17 


SU 


N4450 


Sab 


18.0 


1.2 X 10* 


l068 


gas, 


2 


121±6 


-21.29 






3.46 


SU 


N4477 


SBO:? 


18.0 


8.4 X 10' 


1068 


gas. 


2 


134±6 


-20.92 






1.19 


SU 


N4486B'^ 


El 


17.0 


1.1 X 10' 


l068 


stars. 


6 


185±9 


-17.80 


-17.80±0 


04 


10.1 


SU 


N4501 


Sb 


18.0 


7.9 X 10' 


1CJ68 


gas. 


2 


136±6 


-22.02 






1.51 


SU 


N4548 


SBb 


20.3 


3.4 X 10' 


1068 


gas. 


2 


154±7 


-21.51 






0.71 


SU 


N4698 


Sab 


18.0 


7.6 X 10' 


l068 


gas. 


2 


116±5 


-20.87 






2.13 


SU 


N4800 


Sb 


16.3 


2.1 X 10' 


1068 


gas. 


2 


112±5 








1.00 


SU 


A2052-BCG'' 


E 


151.1 


4.9 X 10' 


logs 


gas. 


7 


233±11<' 


-24.21 


-24.21 ±0 


15 


5.33 


SU 



References. — (1) Pastorini et al. (2007). (2) Sarzi et al. (2002), (3) Gultekin et al. (2009), (4) Marconi et al. (2003), (5) Coccato et al. (2006), (6) Kormendy et al. (1997), (7) 
DallaBonta et al. (2008). 

Note. — is the upper limit to the BH mass at confidence level given by the following column where lOf,^ — 84.2%, 2<7(,s = 97.8%, and 'iO(-,^ = 99.9%. ''Method" column 
indicates method of BH mass determination: stellar dynamics (''stars"), gas dynamics ("gas"). All galaxies in this table are members of the full sample with upper limits (SU). 

^ Bulge-disk decomposition is taken from Lauer et al. (2007b).'^ The galaxy is a pseudobulge.*" NGC 4486B hosts an asymmetric double nucleus (Lauer et al. 1996). Kormendy et al. 

(1997) used ground-based spectroscopic observations of NGC 4486B, along with spherical, isotropic, dynamical models, to find a central dark object of mass between 3 x 10^ and 

3 X lO^M©. Magorrian et al. (1998) used two-integral, axisymmetric dynamical models to find a mass of M = 9.2ljj x lO^M©. The presence of an asymmetric double nucleus coiUd 

indicate an eccentric disk or torus of stars near the center, which requires a massive BH to be stable (Tremaine 1995). The exact mass of the BH , however, may not be sufficiently well 

determined to warrant using the mass found with isotropic models because (1) the three-integral axisymmetric models do not significantly rule out the absence of a BH (Kormendy 
et al. 1997) and (2) the presence of a double nucleus may present problems for isotropic models. To account for this, we conservatively list it as an upper limit at the lOgg upper error 
estimate.*^ Central dispersion Gc is given rather than effective <5e \ see eq. (1). 
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TABLE 7 

Black Hole Masses Omitted from Fits 



Galaxy 




Type 


Dist. 




A/low 


A/high 


Method, 




v,f 


V, bulge 










Mpc 


Mr, 


M^ 


M^ 


Ref. 


km s~^ 






Cygnus 




E 


257.1 


2.7 X 10' 


1.9 X 10' 


3.4 X 10' 


gas, 1 


270 


-21.27 


-21.27 


1.27 


NOIOS* 


MlOl 


Sph 


0.74 






3.8 X lO'* 


stars, 2 


39 


-16.38 


-16.38 


0.03 


N0598'= 


M33 


Sc 


0.80 






3.0 X 10^ 


stars, 3 


24 


-18.77 




0.06 


N3945'' 




SBO+ 


19.9 






5.1 X lO'' 


stars, 4 


192 


-21.06 


-20.09 


1.50 


N4151= 




SAB(rs)ab: 


13.9 


4.5 X lO'' 


4.0 X 10'' 


5.0 X lO'' 


stars, 5 


93 


-20.68 




0.44 


N4303' 


M61 


SABbc 


17.9 


4.5 X lO*" 


2.8 X 10* 


1.4 X lO'' 


gas, 6 


84 


-21.65 




0.31 


N47428 




E4 


16.4 


1.5 X 10' 


9.5 X 10* 


1.9 X lO'' 


stars, 7 


90 


-19.91 


-19.91 


0.99 


N4945'' 




Sc 


3.7 


1.4 X 10^ 


9.0 X 10^ 


2.1 X 10* 


masers, 8 


134 






4.67 


N5252' 




SO 


103.7 


1.0 X 10' 


5.4 X 10* 


2.6 X 10' 


gas, 9 


190 






2.42 



References. — (1) Tadhuntei et al. (2003), (2) Valluri et al. (2005) (3) Gebhardt et al. (2001) and Merritt et al. (2001), (4) Giiltekin et al. (2009), (5) Onken et al. (2007), (6) 
Pastorini et al. (2007), (7) listed a.s in preparation in Tremaine et al. (2002) but never published, (8) Greenhill et al. (1997), (9) Capetti et al. (2005). 

Note. — We list galaxies with dynamically determined BH masses that we do not include in our main fits. We also do not include (but do not list above) the 105 upper limits due 
to Beifiori et al. (2008). These upper limits come from HST spectroscopy of ionized gas at the center of the galaxies, assuming isotropic inclinations of gas disks. If we include all of 
the upper Umits from Beifiori et al. (2008) with our fuU sample, we find a best fit of a = 7.97 ± 0.08, P = 4.41 ± 0.39, Eo = 0.45 ±0.06. In the table footnote for each galaxy we list 
estimates for the parameters of the M-c relation (eq. 2) when that galaxy is included with the rest of ova full sample. For almost all of the exclusions, there is no significant change. 
^ The gas velocity measurement in Cygnus A (Tadhunter et al. 2003) gives strong evidence of the presence of a BH, but the model was not able to explain consistently both their 
infrared Pact and optical [O 111] data, giving no quantitative analysis of the goodness of fit for their models; ot — 8.07 ±0.10, |3 — 4.13 ±0.53, eo — 0.54 ±0.08.^ The upper limit to 
NGC 0205 is considerably below our best-fit ridge line. We omit it from our sample because determination of the mass-to-light ratio at the center of this galaxy may be hampered 
by a very blue nuclear star cluster. It is also a diffuse spheroidal, and we exclude galaxies without either a classical bulge or a pseudobulge. Including this upper limit in our fits 
with the given velocity dispersion of the nucleus, however, does not significantly alter any fit parameters; a = 8. 12 ± 0.08, P = 4.37 ± 0.43, Eo = 0.45 ± 0.06.'^ NGC 0598 is a 
bulgeless disk galaxy (Kormendy & Kennicutt 2004). The M-c and M-L relations apply only to elliptical galaxies or galaxies with a spheroid component. Additionally, it is not clear 
whether to represent the bulge dispersion by the stellar velocity dispersion of the nucleus (as given in this table) or some function of the circular velocity of the disk; ot — 8. 12 ±0.08, 
P = 4.38 ± 0.45, So — 0-45 ± 0.06.^ The mass from NGC 3945, which is a psuedobtilge, is derived from an axisymmetric code (Giiltekin et al. 2009), but the galaxy is a double-barred 
system; a = ±, P = ±, Eo = ±.^ NGC 4151 has been measured by reverberation mapping (Bentz et al. 2006) to have M = 4.57+047 x 10' Mq (after applying a scale factor that 
calibrates reverberation mapping virial products to the M-c relation), in close agreement with the mass derived from the inclined model of Onken et al. (2007) by stellar dynamical 
measurement. The models by Onken et al. (2007) that assume an edge-on inclination, however, cannot rule out M = 0, and the authors cite the noise in their Ax^ contours as cause to 
label it a "tentative estimate." Hence, we do not include this galaxy in our sample; a = 8. 13 ±0.08, P = 4.06±0.40, Eo = 0.46 ±0.06. NGC 4303 has an unresolved point source at 
the center, which is likely an AGN but may also have a significant amount of mass in stars (Pastorini et al. 2007), especially if nuclear stellar clusters are prevalent. The velocity profile 
is also irregular, and the authors do not consider the BH mass estimate completely rehable (Pastorini et al. 2007); a = 8. 13 ±0.08, p = 4.21 ±0.40, Eo = 0.44 ±0.06.^ NGC 4742 
appeared in the sample of Tremaine et al. (2002) hsted as "in preparation," but it has still not appeared in a refereed pubhcation, and thus we do not include it in otxr sample; 
a = 8.13±0.08, p = 4.14±0.39, Eo = 0.44±0.06.'' NGC 4945 has a velocity profile that is asymmetric fl-om one side of the galaxy to the other and no constraint on the disk 
inclination; a = 8. 10 ±0.08, P = 4.36±0.43, Eq = 0.47±0.06.' The best-fit parameter set in the mass modeling produces a reduced of 16.5 (Capetti et al. 2005), indicating apoor 
fit to the data; a = 8. 14 ± 0.08, P = 4.25 ± 0.41, Eo = 0.45 ± 0.06. 



